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1.  INTRODUCTION 


The  Gaowetricai  Dilution  Of  Precision  (GDOP)  concept  is  a 
powerful,  simple  end  widely  used  measure  for  evaluating  the  effective¬ 
ness  of  potential  measurements.  As  its  name  implies,  GDOP  charac¬ 
terises  the  geometrical  effects  present  when  using  multiple  measure¬ 
ments  to  obtain  a  solution  to  the  navigation  problem. 

A  simple  illustration  is  the  determination  of  user  position 
in  a  plane,  from  measurements  of  angles  to  two  kr.ovn  reference 
points.  Position  determination  reduces  to  finding  the  intersection 
of  the  two  straight  lines  emanating  from  the  reference  points  and 
passing  through  the  user's  unknown  position.  If  the  lines  are  nearly 
parallel,  their  intersection  point  is  very  sensitive  to  small 
inaccuracies  in  measuring  angular  directions.  This  is  a  case  of 
poor  GDOP.  More  robust  solutions,  i.e.,  good  GDOP,  results  when 
the  two  lines  are  orthogonal. 

As  important  as  geometry  is,  it  is  not  the  only  consideration. 

It  will  be  shown  that  selecting  measurements  based  solely  on  geometry 
does  not  necessarily  produce  the  optimum  reduction  of  navigation 
uncertainties.  This  study  re-examines  the  question  of  measurement 
selection.  In  addition  to  geometry  the  varying  noise  levels  on 
candidate  measurements  are  considered,  as  well  as  the  a  priori 
navigation  state  uncertainties.  Section  U  considers  generalization 
of  the  GDOP  concept  to  include  these  factors. 

When  a  sequence  of  measurements  are  made  over  a  period  of  time, 
the  optimum  measurements  may  not  be  those  which  give  the  most  accurate 
position  solution  at  each  point  in  time.  This  is  true  for  two  differ¬ 
ent  kinds  of  reasons.  First,  error  growth  between  measurements  is  a 
strong  function  of  velocity  errors,  so  a  position/velocity  error 
trade-off  exists.  Second,  instantaneous  criteria  like  GDOP  do  not 
account  for  visibility  times,  or  lost  measurement  opportunities  dur¬ 
ing  search  and  acquisition  of  new  measurements.  The  result  is  that 
a  saquance  of  less  effective  measurements  over  a  longer  period  of 
time  may  produce  smaller  average  errors  than  would  the  most  effective 
meaaurenents  when  used  over  a  shorter  time. 
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There  is  also  strong  evidence  that  a  law  of  dimenishing 
returns  applies  for  sequential  measurements.  If  the  best  measure¬ 
ment  combination  (baaed  on  Instantaneous  measures  like  GDOP)  is  used 
repeatedly  over  a  sequence  of  measurement  times,  the  incremental 
value  decraasea  with  each  use.  All  other  things  being  equal,  it  is 
possible  that  some  sort  of  altarnation  among  several  combinations 
of  (instantaneously)  less  effective  measurements  will  give  a  smaller 
average  navigation  error.  Section  5  is  devoted  to  some  of  the 
issues  related  to  sequential  measurements. 

If  a  primary  navigation  sensor  fails  completely,  the  system 
is  forced  to  rely  on  some  les9  optimum  mode  of  operation,  perhaps  as 
crude  as  dead  reckoning.  However,  when  a  primary  navigation  system 
is  operational  but  at  a  degraded  accuracy  level,  the  decision  is  not 
so  clear.  Section  6  discusses  some  aspects  of  the  sensor  substitu¬ 
tion  problem.  In  the  Global  Positioning  Satellite  system  (GPS), 
measurements  from  some  satellites  may  be  heavily  jammed,  or  lost 
completely.  The  question  of  sensor  substitution  deals  with  condi¬ 
tions  under  which  a  certain  satellite  channel  ought  to  be  dropped, 
and  if  so,  should  certain  stateB  like  clock  bias  stop  being  updated 
temporarily,  or  should  some  substitute  set  of  inertial  or  radar, 
etc.  hardware  be  used  to  carry  on  in  a  degraded  mode? 

The  scope  of  this  study  contract,  and  the  specific  tasks 
alluded  to  above  are  reiterated  in  Section  3.  Prior  to  this  more 
formal  statement  of  the  problem,  Section  2  provides  some  necessary 
background  development,  including  a  mathematical  review  of  GDOP. 
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2.  BACKGROUND 


This  section  presents  a  brief  tutorial  development  of  GDOP 
from  the  point  of  view  of  least-squares  solutions  to  simultaneous 
linear  equations.  Then  veighted  least-squares  solutions,  with  and 
without  a  priori  statistics,  are  reviewed.  Finally,  the  recursive 
Kalaan  iilter  equations  for  minimum  variant*  estimation  are  presented. 
Tha  well-known  equivalence  (under  certain  conditions)  of  weighted 
least-squares  and  minimum  variance  solutions  is  pointed  out.  Compari¬ 
son  of  these  various  results  will  bring  to  lignt  certain  deficiencies 
in  the  GDOP  approach,  and  will  suggest  generalizations  of  GDOP,  to 
be  developed  in  Section  3. 

2 . 1  An  Introduction  to  GDOP  Using  Least-Squares 

Consider  a  vector  of  measurement  components,  z,  corrupted  by 
additive  noise  n  and  linearly  related  to  the  state  vector  x. 

a  ■  Hx  +  n  (1) 

An  estimate  of  x,  called  x,  is  sought  from  the  given  z  vector.  One 
approach  to  this  problem  is  the  least-squares  estimator  x.  T^at  is, 
select  x  ao  as  to  minimize  the  sum  of  the  squares  of  the  components 

A  * 

of  the  difference  z  -  Hx.  This  means  x  should  minimize 

Jj  ■  (z  -  Hx)T(z  ~  Hx)  (2) 

Differentiating  with  respect  to  x  and  setting  the  result  to  zero 
gives 

T  T  * 

-  H  z  +  H  Hx  -  0 

so  that  the  well-known  least-squares  estimate  is 

x  -  (HTH) -1  HTz  (3) 


The  quality  of  this  estimate  can  be  assessed  by  forming  the  error 

x  -  x  -  x  (A) 


Using  Eq.  (1)  and  Eq.  (2)  gives 

X  -  (HTH)_1  HTHx  +  (HTH)-1  ’lTn  -  x 

-  (HTH)_1  HTn  (3) 
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If  the  noise  n  has  itro  mean,  taking  the  expected  value  £( }  of 
Sq.  (5),  shows  that  the  estimator  error  x  also  has  a  zero  mean. 

The  covariance  of  the  estimation  error  can  be  related  to  the  covari¬ 
ance  of  the  noise  n. 

E(x  xT>  -  (HTH)-1  HT  E{ n  nT}H(HTH)-i  (6) 


If  all  components  of  n  are  pairwise  uncorrelated  and  have  unit 
variance,  i.e., 


then 
so  that 


*Hni  n^} 

E{n  nT} 


if  i  -  j 
if  i  +  j 


1 


E{x  x  } 


T  -1  T  T  -1 
(HH)  H  H(H  H) 

T  -1 
(HH) 


(/) 


T  -1 

The  matrix  (H  H)  will  be  called  the  GDOP  matrix.  Various  scalar 

terms  derived  from  this  matrix  are  frequently  used.  For  example, 

T  -1  t  -1 

the  trace  of  (H  H) ,  written  TR(H  H) ,  is  the  sum  of  the  squares 

of  all  error  components  when  E{n^  n^ }  -  5  .  The  square  root, 

[TR(HTH)  1]ls  is  usually  called  GDOP.  ^ 1 ^  The  square  root  of  the 

sum  of  all  diagonal  terms  which  correspond  to  position  components 

r  i  2 1 

is  sometimes  called  PDOP. 1  ’  J  If  the  horizontal  plane  diagonal 

t  *1  r i  2 1 

components  of  (H  H)  are  summed,  the  square  root  is  called  HDOP.  ’  1 

T  -1 

In  this  study,  major  attention  is  directed  to  TR(H  H)  .  Omission  of 
the  aq tare  root  will  not  alter  any  conclusions  regarding  the  minimiza¬ 
tion  of  GDOP. 

It  is  clear  that  all  GDOP-related  performance  measures  indicate 
the  error  in  an  estimated  navigation  quantity  "per  unit  of  measurement 
noiBe"  covariance.  Deviations  from  this  assumption  are  considered  in 
the  next  Section. 

All  of  the  above  GDOP-related  measures  depend  solely  on  the 
geometry  matrix  H.  Smaller  GDOP  values  indicate  stronger  or  more 
robust  geometric  solutions  to  the  estimation  problem.  For  these 
reasons,  when  some  freedom  exists  in  the  choice  of  measurements, 
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good  (i.e.,  small)  GDOP  Is  often  used  as  tha  selection  criterion. 

It  vlll  be  ahovn  later  that  this  Is  not  alvays  the  best  choice. 

2 . 2  Weighted  Least-Squares  Solution 

The  "per  unit  noise"  concept  inherent  in  GDOP  is  not  so  useful 

(or  more  bluntly,  invalid)  when  certain  measurement  components  are 

noieier  than  others.  When  a  choice  exists  betveen  two  possible  sets 

of  measurements,  the  set  with  poorer  GDOP  may  be  preferable  if  they 

are  of  sufficiently  higher  accuracy.  This  trade-off  was  pointed  out 

U  | 

in  conjunction  with  an  integrated  GPS/Inertiai  Navigation  System. 1,}J 
There,  the  desire  to  achieve  good  GDOP  suggested  selection  of  range 
measurements  separated  by  angles  near  90°.  However,  this  caused  an 
increase  in  atmospheric  diffraction  errors  in  those  measurements 
with  low  elevation  angles.  Another  potentially  more  severe  problem 
with  nonuniform  noise  arises  in  a  jamming  environment.  Some  of  the 
highly  directional  antennas  may  receive  large  amounts  of  noise 
pollution.  This  degradation  would  more  than  overcome  any  GDOP 
advantages . 

If  the  measurement  noise  covariance  is  not  just  the  unit 

T  2 

matrix,  but  rather  is  E{n  n  }  •  Io  then  the  GDOP  matrix  just  ratios 
2 

up  by  the  scalar  o  .  If  all  potential  measurements  have  this  same 
2 

variance  o  ,  the  choice  of  the  best  measurement  set  will  still  be 

2 

the  best  GDOP  set.  The  multiplication  scalar  a  will  not  affect 
the  relative  rankings, 

Tha  measurement  noise  is  called  nonunifonn  when  different 

measurements  have  different  noise  levels,  as  indicated  by  their 

U1 

variance*.  In  this  case  a  weighted  least-squares  approach1-  J  to 
estimation  is  often  used.  The  quadratic  form  of  Eq.  (2)  is  modified 
by  inserting  a  weighting  matrix  W, 

J2  -  (i  -  Hx)T  W(r  -  Hx)  (8) 

It  is  customary  to  select  the  weighting  matrix  as  the  inverse  of  the 

-1  T 

noise  covariance  matrix  (i,e..  W  •  R  where  R  ■  E{n  n  }).  This 
choice  weights  tha  accurate  measurements  more,  the  noisy  ones  less. 
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Minimizing  this  modified  cost  function  lends  to  a  modification  of 
Eq.  (2) 

*  T  -1  -1  T  -1 

x  -  (H  R  H)  H  R  z  (9) 

and  a  corresponding  .uodif ication  of  £q.  (6)  or  (7) 

Elx  xT}  -  (HTr"1H)~1  (10) 

Equation  (10),  and  generalizations  of  it,  will  be  discussed  in  termB 
of  potential  CDOP  generalizations  in  Section  A.  It  is  presented 
here  only  to  contrast  the  traditional  GDOP  with  the  resulting 
estimation  *rror  covariance. 

Thus  far  no  mention  has  been  made  of  the  knowledge  about  the 
navigation  states  which  may  exist  prior  to  making  navigation 
measurements.  In  almost  all  real  situations  prior  knowledge  will 
exist  because  of  previous  measurements  or  for  other  reasons.  The 
uncertainty  volume  for  x  need  nor.  be  sphericul.  Some  components  of 
x  esn  have  much  more  uncertainty  associated  with  them,  leading  to 
elongated  ellipsoidal  uncertainty  volumes.  Those  components  with 
the  largest  uncertainty  are  said  to  have  "maximum  need"  for  improve¬ 
ment.  Some  measurements  with  less  than  optimal  over-all  GDOP,  may 
be  more  effective  in  reducing  navigation  uncertainties. 

The  treatment  of  a  prior  information  in  estimation  cun  be 
approached  in  several  ways.  One  way1  J  is  to  again  modify  Eq.  (2), 
this  time  to 

•  (z  -  Hx)^  R  ^"(z  -  iix)  +  (x^-x)'*'  ?0^(xo_  *)  (ID 

where  x  represents  an  a  priori  estimate  of  x  and  the  covariance  of 

— o  — 

x  Is  P  .  Minimizing  J*  amounts  to  a  trade-off  between  measurement 
-00  3 

residuals  and  deviations  from  x  .  This  leads  to  a  modified  estimator 

-o 

x  -  (HTR_1H  +  P_1]”1tHTR“1z  +  p”1  x  ]  (12s) 

_  k  O  —  o  -o 

-  x  +  P  HT[HP  HT  +  R]-1[z  -  Hx  ]  (12b) 

— o  0  0  —  o 

Form  (12a)  most  commonly  is  U3ed  in  weignted  least-squares  analysis, 
while  for*  (12b)  is  ons  of  the  well-known  Kalman  filter  equations^  . 
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These  two  form*  are  in -fact  identities.  ,  and  tha  corresponding 
forms  of  ths  a  posts- iori  estimation  error  covariance  matrix  are 
givan  by  the  following  identities . ^ 

E{x  xT}  -  PJ  -  tHTR_1il  +  p"1]"1  (13a) 

-  Po  -  PoHT[H?0HT  +  R]"1  HPo  (13b) 

Mote  that  Eq.  (13a)  loots  most  like  the  GDOP  matrix  in  that  it  reduces 
to  F.q.  (10)  If.  (very  large  initial  uncertainties)  and  it 
reduces  to  Eq.  (7)  if,  in  addition,  R  *  I.  These  equations  will  be 
used  in  Section  3,  when  discussing  a  generalised  GDOP  concept. 

They  are  presented  here  to  enable  the  deficiencies  of  the  traditional 
GDOP  sethod  to  be  demonstrated. 

2.3  An  introductory  Example 

Consider  a  planar  navigation  problem.  Range  measurements  from 
the  unknown  vehicle  location  x  to  two  known  points  x^  and  x^  ere  used. 
The  geometry  Is  shown  in  Figure  I.  The  range  measurements  are  and  z. 


Figura  1.  Gsomstry  of  a  Flsnsr  Rsngs  Navigation  System 


are 


Neglecting  measurement  noiea,  tha  rang** 


.«•  /  <*-*,> 2  * 


)■ 


(14) 


S  "  /  +  <y'yb)2 

These  for*  components  of  the  nonlinear  measurement  aquation 

z  •  h(x) 

In  ordar  to  uaa  <3)0?  in  thia  exerple,  linaar  sensitivity  aquations 
ara  required.  Linearizing  about  some  nominal  point  xq  gives 


Az 


3t  3z 

Ti  *x  +  T?  iy 


a*. 


I  W4,K 

I  Ax  +  ~  Ay 

L  3x  ay 


V*.  Vya 


a  a 

vs  vs 


Ax 


L  Ay  J 


or  A*_  -  H  Ax 
The  component*  in  H  c 
geometry  of  Figure  1, 


(16) 

(17) 


ba  written 

as  trig 

sin  9 

cos  6 

a 

* 

sin  9b 

cos  8b 

(18) 


The  GDOP  matrix  i*  expressible  us 

(hth)'l  “ 


T 


■in  (ea-«b) 


2  2 
coa  6^  +  cos 


4(sin  29^+  sin  26b) 

2.  .  .1 


-Jj(sin  26a  +  sin  26^  *in  9^  +  sin  6fe 
Note  that  the  planar  GDOP  is 


(19) 


■nto^H)"1 


(20) 


sin  (0-9. ) 

•  D 
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Tbit  la  independent  of  th«  choice  of  coordinate  systems  and  is 


min  learn  (beat)  when  e  -  6.  - 
the  off-diagonal  terms  in  (H 


±  90°.  When  0  -  8,  •  t  90°,  note  that 

T  -1  *  ® 

H)  are  aero,  indicating  a  lack  of 


cross  coupling  of  errors •  The  geometrical  interpretation  of  these 
results  ie  clear.  The  point  x  -  is  being  determined  by  the 
intersection  of  two  circles  of  radii  z^  and  z^" respectively, 
centered  at  x#  and  x^.  Small  errors  in  measuring  or  z^  have 
the  least  effect  when  the  circles  intersect  at  right  angles. 


-If-0^--  8^  is  small,  so  the.  circles  graze  each  other  at  a  shallow 
angle,  then  smell  errors  in  z  or  z  give  large  errors  in  x.  An 

H  D  ** 

infinite  error  results  if  6  ■  0,  .  This  means  there  is  no  inter- 

a  b 

section  at  all  if  t  it.  and  no  unique  intersection  if  z  ■  z.  . 

AD  AD 

Soma  spaclfic  numerical  cases  related  to  Figure  1  are 


informative.  Suppose  there  are  two  candidate  sets  of  measurements 


Sat  1*6  •  90°,  0.  -  0  ,  U 

a  b 


Sat  2i  0  -  63.435,0  -  -45°  ,  H 

a  b 


From  Eq.  (20)  tha  trace  of  the  GDOP  matrix  for  Sec  1  is  2.00. 

This  ie  the  minimum  possible  value  and  5et  1  is  therefore  optimal 
in  terms  of  GDOP .  Sec  2  has  a  corrasponding  suboptlmal  value  of 
2.222.  Tht  optimality  of  Set  1  is  expected  since  0  -  6^  -  -90°. 

For  Set  2,  0.-  0  •  -108.435°  so  the  unit  measurement  vectors  are 

D  A 

not  quite  orthogonal,  giving  the  suboptimal  GDOP. 

Suppose  that  these  two  sets  of  measurements  are  to  be  compared 
in  terms  of  tha  noise  weighted  GDOP  matrix  of  Eq.  (10).  For  Sat  1, 
•inca  H  -  I, 

E{x  xT)  -  R 
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Suppose  chat  Sat  1  has  R  ■ 


4.5 

0 


0 

2.5 


and 


Sat  2  has  R 


4  0 

,0  2 

Than  the  trace  ol  Eq.  (10)  for  Sat  1  has  a  value  of  7.  For  Sat  2, 
tha  corresponding  value  is  20/3  -  6.667.  A  vary  modest  decrease 
in  noise  variances  was  sufficient  to  overcome  the  GDOP  advantage 
and  make  Set  2  the  best  choice  in  terms  of  the  weighted  least- 
squares  ~r iter ion. 

In  order  to  demonstrate  the  effect  of  a  priori  information, 

assume  that 

'3  -1 

P  - 

°  -1  3 


and  compare  measurements  Sets  1  and  2  in  terms  of  the  trace  of  the 
a  posteriori  covariance  in  Eq.  (13).  To  nullify  effects  of  different 

^measurement  noise  levels,  assume  that 


R  - 


f 

4 

.0 


0 

2  . 


for  both  measurement  sets.  Straightforward  calculation  of  Eq.  (13a) 
shews  that  T)«x  xT)]  is  2.823  for  Set  1  and  2.697  for  Set  2. 

Thus  Set  2  gives  the  smaller  a  posteriori  error  volume  even  though 
its  GDOP  value  is  inferior  to  Set  1  and  even  though  the  same  measure¬ 
ment  accuracies  are  assumed  in  both  cases.  These  simple  results 
suggest  that  minimum  GDOP  does  not  guarantee  minimum  a  posteriori 
position  uncertainties  in  many  cases.  Similar  results  are  found 
in  realistic  GPB  examples  to  be  presented  later. 


2.4  Kiniwum  Variance  Estimation 


If  instead  of  seeking  to  minimize  a  selected  cost  function 
such  as  found  in  Eq.  (2),  (8)  or  (11),  one  seeks  the  estimate  x 

T  A 

which  minimizes  the  error  variance  E{ (x-x)  (x-x) } ,  then  the  minimum 
variance  estimator  will  be  obtained.  Under  quite  general  conditions 
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the  solution  is  ths  conditional  mean 

x  -  E{x| z)  (21) 

If  it  is  a»*u**d  that  the  state  x  and  tha  measurement  x  ar«  jointly 
distributed  gaussian  random  variables,  or  if  the  best  linear 
estimator  is  sought  regardless  of  the  form  of  the  distribution 
functions,  than  Eq.  (21)  is  known  to  yield 

x  -  E{x)  +  ?„  P"J[i  -  E{  z }]  (22) 

If  the  linear  measurements  of  Eq.  (1)  are  assumed,  it  is  well  known 
that  Eq.  (22)  reduces  to  Eq.  (12b),  with  xq  ■  E{x),  The  corresponding 
a  posteriori  error  covariance  matrix  is  given  by  Eq.  (13b),  so  that 
the  weighted  least-squares  and  minimum  variance  solutions  are  the 
same  under  appropriate  conditions. 

It  will  be  useful  for  future  discussions  on  time-sequential 
measurements  to  list  the  other  well-known  equations  which  make  up 
the  recursive  form  of  the  minimum  variance  estimator  known  as  the 
Kalman  filter. ^  The  notation  used  earlier,  such  as  in  Eq.  (13), 
with  subscript  zero  indicating  conditions  Just  prior  to  a  measure¬ 
ment  and  subscript  one  indicating  conditions  Just  after,  is 
inadequate  whan  time-sequences  are  considered.  Two  time  Indices  are 
useful.  The  first  index  indicates  the  time  at  which  the  estimate 
applies  and  the  other  indicates  the  most  recent  measurement  used  in 
forming  the  estimate.  Equation  (12b)  becomes,  in  this  standard 
Kalman  filter  notation 

x (k+l| k+1)  -  x(k+l | k)  +  P(k+l|k)H(k+l)T[H(k+l)P(k+l|k)HT(k+i) 

+R(k+l)]"1[z(k+l)  -  H(k+l)x(k+l|  k) ]  (23) 

Tne  covariance  update  equation  (13b)  becomes 
P (k+1  k+1)  -  P (k+lj  k) 

-P(k+l|  k)HT(k+l) [H(k+l)P(k+l j k)HT+R(k+l) ]_1H(k+l)P (k+l|  k) 

(24) 

The  two  new  equation#  give  the  extrapolation  between  measurement 
times,  and  are 

x(k+l|k)  -  ♦(k+l,k)x(k|k)  +  T(k+1 ,k)E{w(k) }  (25) 
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(26) 


and  P(k+ljk)  -  ♦(k+l,k)P(k|k)eT<k+l,k)  +  Q(k) 

wh«t*  ♦(k+l,k)  ie  ths  state  transition  matrix 

B{v(k)}  la  tht  »aan  plant  noise,  often  aero 
r(k+-l,k)  is  the  input  coupling  matrix .  _ 

*ad  Q(k)  is  the  covariance  matrix  of  the  plant  noise 

r(k+l,k)w(k)  -  . 


■)  ■ .  urn  wrmOL ** 


3 .  SCOPE  OP  THE  STUDY 


Tbit  cuuich  study  addresses  three  major  tasks,  S3 

outlined  belov. 

3.1  Measurement  Selection  Methodology 

Based  upon  the  results  of  the  simple  example  in  Section  2.3, 
and  previous  experience,  It  is  known  that  the  traditional  GDOP 
does  not  always  indicate  which  sets  of  potential  measurements  are 
optimal.  In  addition  to  the  geometry  effects,  which  are  accounted 
for  by  GDOP,  nonuniform  measurement  noise  and  a  priori  state 
information  is  taken  into  account  and  a  scalar  GDOP-like  indicator 
of  static  measurement  effectiveness  is  presented. 

When  a  sequence  of  measurements  are  to  be  taken  over  a  period 
of  time,  the  dynamics  of  the  system  and  the  operational  scenarios 
become  factora  in  measurement  selection.  This  aspect  of  measure¬ 
ment  selection  methodology  la  also  addressed. 

If  a  portion  of  the  primary  navigation  system  1b  degraded 
because  of  jamming  or  other  reasons,  then  secondary  or  alternate 
sensors  must  be  considered.  This  aspect  of  the  measurement  selection 
strategy  is  addressed,  using  auboptimal  filtering  and  partitioned 
state  vector  cone  epts. 

3.2  Analytical  Tests 

The  methods  which  result  from  the  task  of  Section  3.1  are 
evaluated  to  the  extent  practical  by  analytical  methods.  Simple 
models  are  used  for  gaining  insight,  and  realistic  models  are  also 
analyzed.  Numerical  results  are  also  presented  for  realistic  GPS 
geometries  in  some  cases.  It  is  pointed  out  the  large  scale  simula¬ 
tions  were  explicitly  excluded  from  this  study  because  of  the  costs 
that  would  have  been  involved  in  developing  the  necessary  capability 
which  is  alrsady  in  existence  at  the  Air  Force  Avionics  Laboratory. 
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r 


3.3  Specification  of  Simulation  Casas  and  Condition  a 


.  A*  mentioned  In  Sacclon  3.2,  the  Intent  of  this  study  was 
not  to  duplicate  existing  siaulation  capabilities.  This  study 
does  suggest  certain  types  of  cases  and  conditions  that  ought  to 
be  considered,  via  siaulation,  in  further  evaluation  of  the  concepts 
developed  here. 
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4.  GENERALIZATIONS  AND  EXTENSIONS  OP  GDOP 


There  1*  absolutely  no  quticlon  chat  tha  GDOP  concept  has 
playad,  and  will  continue  to  play  a  major  role  in  the  design  of 
many  navigation  systems.  It  is  simple  to  use  and  haB  modest 
computational  requirements.  It  cheraccerltes  tha  geometrical 
aspects  of  a  navigation  satellite  network  very  nicely.  In  the 
case  of  satellite  network  design,  the  basic  requirement  is  to 
assure  that  a  vide  variety  of  geographically  dispersed  users, 
with  a  vide  variety  of  sensor  accuracies  and  a  priori  knowledge 
are  not  prevented  from  solving  their  navigation  problem  because 
of  veak  geometry. 

It  is  when  GDOP  is  proposed  as  a  means  for  an  individual 
user  to  select  measurements  which  are  best  under  his  particular 
circumstances  that  room  exists  for  improvements,  generalizations 
and  extensions.  A  major  part  of  this  study  was  devoted  to  the 
question  of  how  to  best  select  measurements  at  a  given  point  in 
time.  This  statie  measurement  selection  problem,  so  called 
because  only  one  point  in  time  is  considered,  is  discussed  in 
this  section. 

To  begin,  a  definition  of  what  is  meant  by  "best"  is  desirable. 
The  point  of  view  is  assumed  here  that  best  means  the  smallest  over¬ 
all  mean  square  position  error.  In  the  GPS  problem  we  are  concerned 
with  three  position  components  and  a  clock  bias,  which  converts  to 
en  equivalent  range  error.  So  for  the  most  part  a  four-component 
state  vector  will  be  of  ultimate  interest.  Making  the  usual 
assumptions  about  mean  values  being  zero,  attention  turns  to  the 
second  statistical  moments  of  the  estimation  errors,  l.e.,  the 
covariance  terms.  Subject  to  the  uausi  assumptions,  Eq .  (13a)  or 
(13b)  or  Eq.  (24)  tells  the  entire  story.  Given  an  a  priori  state 
uncertainty  Pq,  and  a  sat  of  measurements  characterized  by  a 
geometry  matrix  H  and  a  noise  covariance  matrix  R,  the  a  posteriori 
covariance  is  uniquely  computable. 
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Ic  ia  awkward  to  uaa  a  A  x  4  matrix  to  describe  the  accuracy 
or  ef factivantas  of  a  given  eat  of  measurements.  An  attractive 
alternative  la  to  uaa  a  single  scalar  number  to  deacrlbe  the  result. 
In  the  past,  various  terns  such  as  circular  error  probability  (CEP) 
or  spherical  error  probability  (SEP)  have  been  used. 

The  a  posteriori  covariance  matrix  P  can  be  viewed  geo¬ 
metrically  as  a  hyperellipsoid  with  principal  axes  oriented  in 
space  according  to  tha  diractions  of  the  orthonoraal  eigenvectors 
of  P^  and  with  the  dimensions  of  the  semi-major  axes  as  the 
eigenvalues  X^  of  P^.  The  interior  of  this  ellipsoid  can  be  thought 
of  as  a  ona-sigma  error  volume,  and  in  this  respect  it  is  related 
to  CEP  and  SEP  although  much  more  complicated.  Two  very  simple 
scalar  maasuras  derived  from  P1  have  been  used.  The  determinant 
of  P^,  is  equal  to  the  product  of  the  eigenvalues^  J 


P  m  \ 

1  A1 


X «  ...  X 

2  n 


and  is  therefora  ralatad  to  the  volume  of  the  ellipsoid.  This 
measure  is  not  preferred  here  for  two  reasons.  First,  determinants 
are  not  as  easy  to  evaluate  as  the  trace,  which  is  used  in  the 
next  performance  measure.  Second,  and  more  important,  | |-  0  does 
not  mean  the  error  is  zero.  All  it  means  1b  that  one  dimension  of 
the  hyperelliosoid  has  degenerated  to  zero,  but  very  large  errors 
in  the  other  coordinates  may  still  exist.  |  P^  | is  not  very  descriptive 
or  discriminating  in  this  regard. 

The  preferred  scalar  performance  measure  is  based  on  the  trace 


TR(P  )  -  Ip}1  -  l  X 
i-l  i-l 

The  trace  will  normally  ba  computed  as  the  sum  of  the  diagonal 
elements,  but  it  is  known  that  it  also  equals  the  sum  of  the 
eigenvalue* . ^  Therefore  TR(P^)  can'not  equal  zero  unless  all 
error  variance  components  era  zero. 
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There  may  be  circvautancee  where  it  it  sort  Appropriate  to 
uee  e  veighted  eua  of  the  diagonal  terms  of  P^  as  a  performance 
. measure.  The  vertical,  component  of  position  error  nay  be  more 
important  than  the  horirontal  components  In  a  terrain-following 
ayatea,  for  example.  By  defining  a  diagonal  matrix, 


(29) 


It  Is  easily  verified  that 

n  n 

TR(WP.)  -  l  a  P  (30) 

i-1 

For  many  purposes  TR(WP^)  can  be  dealt  with  just  as  easily  as  TR(P^) , 
and  they  are  the  same  when  W  is  the  unit  matrix. 

The  trace  functions  just  defined  are  natural  extensions  of 
GDOP.  In  the  remainder  of  this  section  several  versions  of  extended 
GDOP  will  be  considered,  Analytical  techniques  using  eigenvalue- 
eigenvector  theory  are  developed  in  an  attempt  to  gain  insight  Into 
optimal  measurement  selection.  Computational  considerations  prompt 
the  evaluation  of  some  approximate  but  easy  to  compute  measures. 
Numerical  rankings  of  realistic  measurement  possibilities  are 
provided,  to  show  the  performance  of  the  various  exact  and  approximate 
methods . 


4.1  Noise-Weighted  GDOP 

In  Section  2.1  it  wee  agreed  that  the  trace  of  Eq.  (7)  would 
be  called  GDOP.  The  trace  of  Eq.  (10)  is  called  the  noise-weighted 
GDOP  for  obvious  reasons, 

Noiae-veighted  GDOP  -  TR(HIR~1B)''1  (31) 

The  Intent  is  to  use  measurement  accuracy  information  as  veil  as 
geometry.  This  performance  measure  will  prevent  the  selection  of 
very  noisy  (large  R)  measurement  sets.  This  rather  minor  modifica¬ 
tion  to  GDOP  is  often  effective.  In  feet,  comparison  with  Eq.  (13a) 
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shows  It  to  be  exactly  TR(P.)  for  the  case  where  uncertainties  in 

-1  1 

PQ  ere  to  lerge  that  Pq  can  be  neglected.  While  this  circumstance 
»ay  apply  in  the  early  stages  of  a  navigation  solution,  it  is  not 
particularly  appropriate  for  a  well-behaved  navigation  system  in 
steady-state  operation.  It  will  be  shown  in  numerical  examples  later 
that  noige-veighted  GDOP  can  occasionally  do  worse  than  traditional 
GDOP  in  selection  of  measurements.  Mora  often  it  is  better  than 
GDOP  but  still  suboptimal. 

4.2  Incorporation  of  A  Priori  Dsta 

Both  Eq.  (13a)  and  (13b)  ve  the  a  posteriori  covariance. 

Eq.  (13a)  is  most  similar  to  thu  noise-weighted  GDOP  matrix,  differing 

only  by  the  inclusion  of  the  a  priori  covariance  P  .  However, 

o 

Eq.  (13b)  seems  most  convenient  and  will  be  used  here. 

Consider  as  the  basic  indicator  of  measurement  effectiveness 

TRCWP. )  -  TRCWF  )  -  TR{ WP  HT[HP  HT+  R]_1HP  }  (32) 

X  O  0  0  0 

Since  W  and  are  unaffected  by  choice  of  measurement,  Eq.  (32)  can 
ba  minimized  (with  respect  to  the  choice  of  H  and  R)  by  maximizing 

SCORE  -  TRtVP  HT[H?  HT+  Rj"1  HP  }  (33) 

o  1  o  J  o 

-  TR{HPoWPoHT[HPoHT+  R]"1}  (34) 

Equation  (34)  follows  from  Eq.  (33)  by  using  properties  of  th>  trace. 

The  score  is  the  amount  of  uncertainty  reduction  achieved  from  a 
given  measurement  set.  If  the  score  is  computed  for  each  possible 
measurement  set,  the  maximum  score  will  indicate  the  optimum  selection, 
When  four  satellites  are  to  be  selected  out  of  seven  assumed  to  be 
available  for  measurement,  35  combinations  must  be  evaluated.  In 
earlier  studies,  up  to  ten  satellites  were  assumed  visible,  giving 

210  combinations  of  four.  These  are  large  numbers  of  possibilities 

to  compute  and  rank.  Therefore  effective  algorithms  for,  or  good 
approximations  to,  Eq.  (34)  are  essential.  These  algorithmic  aspects 
will  ba  pursued  in  Section  4.4.  First,  in  order  to  gain  insight  and 


aid  in  developing  *£11018111  ruleB  of  thumb,  an  analytical  Investi¬ 
gation  la  carried  out  using  eigenvalue-  igenveetor  theory. 

4.3  Eigenvalue-Eigenvector  Analysis^ 

The  n  x  n  real  symmetric  a  priori  covariance  matrix  has  a 
full  aet  of  n  real  orthonormal  eigenvectors  5  and  n  corresponding 
real  nonnegative  eigenvalues.  They  saticfv  the  eigenvalue  equation 


Po  Si."  ^  ^ 


i  *i 


i  ■  1 , 2 , . . .  ,n 


(35) 


The  computation  of  eigenvalues  and  eigenvectors  for  a  large  matrix 
can  be  a  burden,  and  is  best  avoided  in  any  real-time  scheme  such 
as  measurement  selection.  However  they  provide  answers  to  interesting 
theoretical  questions  and  are  therefore  considered  in  this  Section. 

The  matrix  Po  can  be  written  in  the  so-called  spectral 
decomposition  form 


P 

o 


n 

y 


ittl 


X1  h  si 


(36) 


4.3.1  Individual  Scalar  Measurements. 

The  effect  of  a  single  scalar  measurement  is  considered  first. 


z  ■  h  x  +  n  (37) 

This  is  a  simplified  version  of  Eq.  (1)  but  now  z  and  n  are  scalars 
and  h  is  a  1  x  n  row  vector.  When  Eq.  (37)  is  used,  the  simplified 
version  of  Eq.  (13b)  is 

P  -  P  -  P  hTh/(h  P  hT  4  R)  (38) 

1  0  Q---0- 

The  measurement  noise  n  has  the  scalar  covariance  R. 

The  set  of  eigenvectors  (5.,  i  ■  l,...,n}  can  be  used  as  a 

T 

basis  for  an  n-dinensional  vector  space,  so  the  column  vector  h  can 

be  expanded  as 

T  n 

h  -  [a  5  (39) 

1-1 

where  the  expansion  coefficients  are  given  by  the  scalar  product 

ai  -  b  M  Si  1)T  (40) 
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Frasultiplyiag  <  (39)  by  PQ  a ad  using  Eq.  (35)  gives 


-  -  ,  •  Pe  *  “  J,X  1  ai  V 

Multiplication  of  Eq.  (41)  by  h  and  using  Eq.  (40)  givas  the 
denominator  required  in  Eq.  (38). 

£  V&T  +  R  "  I  xi  3i  +  R 

i-1  1  --- 

Tha  spectral  fora  of  Eq.  (38)  Is  found  to  be 


(41) 


(42) 


T  <i\”i  $l)(S‘kW 


.  i  h  i.  si  - " '  '  '  c;  ~ 

S  *3  + 11 

By  separating  out  thoie  term*  in  the  double  summation  for  which 
i"k,  this  bacoaes 

°  '  Xi  “i  ->  T 

Pl- 

i-ll  5  Xj  aj  +  R  J 


i^lp,  ai  +  R' 


(44) 


i*k 

i-k 


J  J 

Ortb  mormality  of  the  eigenvectors  gives 

T  t  f  ® 

ra<5i  Sk)  -  ™‘£k  £i>  ■  \i 

Tharafora  the  trace  of  Eq.  (44)  is 

IMP )  -  I  {»  -  —^5 -  \ 

a  5  x3  +  *  > 

The  score  defined  in  Eq.  (33)  and  (34)  is  found  from  Eq.  (45)  for  a 
scalar  measurement  and  with  W  -  I, 


(45) 


2  2 
n  a  a . 

SCORE  -  l  - - - =j - 

i-1  (J  Xj  +  R) 


(46) 
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When  complete  freedom  of  choice  ie  eveileble  regarding  the  direction 

of  the  measurement  vector,  It  la  found  that  Eq.  (46)  la  maximized 

T 

-  4d  Sq.  (45)  »inimi«ed  if  h  -  the  eigenvector  with  Che  largest 
eigenvalue. 

\2 


■in  TR(P, ) 
h  1 


i*<p0)  -  r 


■ax 


+  R 


(47) 


When  the  noiae  covariance  R,  la  allowed  to  vary,  but  the  directions 
T  1 

of  are  reatricted  to  eigenvectora,  the  optimal  measurement  is 
atlll  the  direction  of  maximum  eigenvector  if 


Ri  ^ 


(X.-  X  )  +  R  (X1/X  ) 

— rbrf — 


for  each  other  candidate 


j 

X^  and  . 


Clearly,  if  Xi  ■ 


(48) 


Xj  this  indicates 


In  general  a  leas  accurate  measure* 


that  the  smallest  la  optimal. 

ment  (larger  R^)  may  be  optimal  if  there  is  more  to  be  gained  (X^  s  x  ) 

Whenever  a  scalar  measurement  vector  h  is  in  the  direction 

of  any  eigenvector,  then  there  is  only  one  nonzero  a  expansion 

T  1 

coefficient  in  Eq.  (39).  Assume  that  h  Is  a  unit  vector.  Then 
<*r  »  1  ,  •  0  j^r  and  Eq.  (44)  becomes 

2  2 


1 '  Ji1'  ‘  fers]  5“)  Sl 

'  ill  «+  h  hr’}  S‘  ^ 


(49) 


where  6 


ir 


The  implications  of  this  result  are  that  P^  has  the  same 

eigenvectors  as  Pq  and  all  the  eigenvalues  ars  unchanged  except  Xf, 

In  terms  of  the  error  hyperellipsoid,  the  orientation  of  its  axes 

remains  fixed  in  apace.  All  coordinate  dimensions  sre  also  the  same 

as  before  the  measurement,  except  the  one  in  the  direction  of  the 

measurement  is  squeezed  down  from  X  to  X  R/(X  +  R) .  If  R  -  0,  this 

r  r  r 

dimension  flattens  to  zero.  If  R  Is  larg<>,  virtually  no  reduction 
takes  place. 
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A  aore  likely  decision  to  be  faced  la  the  selection  of  the 
best  aeasureaent  from  aaong  several  vhich  are  not  necessarily  in 
the  direction  of  eigenvectors  of  Pq.  A  straightforward  numerical 
decision  based  on  Eq.  (34)  will  be  presented  In  Sections  4.4  and 
4.3.  Here,  a  tvo-diaenslona]  case  is  considered  analytically. 

As suae  that  the  eigenvectors  and  eigenvalues  define  the 
error  ellipse  shown  in  figure  2. 


Figure  2.  Error  Ellipse  and  Two  Candidate 
Measureaent  Directions 

The  candidate  aaasureaent  vectors  make  angles  &  and  y  respectively 
with  the  sealaejor  axis.  Noise  covariances  end  R^  *r*  assumed. 
The  ecelar  product  of  two  unit  vectors  is  tbs  cosios  of  thslr 
lncludsd  sngls, 

hx  •  cos  6,  h2  ^  -  cos  y 

Thsrsfors  Eq.  (46)  can  bs  expressed  as 
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(50) 


and 


2  2  2  2 
X.  cob  g  +  A?  sin  B 

3  ■  .  r  - .  .  .. 

*  A.  coa^B  +  A 


S2  " 


, 2  2  ,2  2 

A^  cos  y  +  Aj  *in  y 

2  2 

A2  cob  y  +  *2  #in  Y  +  R2 


for  tha  tvo  potential  auBurnwnti. 

In  order  co  determine  whether  or  S2  ie  larger,  tha  difference 
is  computed,  a  positive  common  denominator  is  omitted  and  the 
results  manipulated  into 


A  positive  d  indicates  that  choice  1  is  best.  Several  special 
cases  are  of  interest.  If  A2/A^  -  1  so  the  ellipse  is  actually  a 
circle,  than  tha  measurement  with  the  smallest  R1  is  best  regardless 
of  the  angles  y  and  6.  If  A2/A^  **  negligible,  so  we  have  an 
extremely  eccentric  ellipse,  it  can  be  shown  that  choice  1  is  best  if 

2  2 

cos  y  <  R2  C0B  8  (51) 

One  ocher  special  case  illustrates  the  trade-offs  between  eccentricity, 

T 

off-axis  angles  and  noise  levels.  Let  y  ■  0  so  chat  Is  along 

the  maximum  eigenvector.  It  can  be  shown  chat  d  remains  positive, 

X 

an  that  the  off-axis  measurement  h^  is  still  best  if 


/  (R,-  R  )/R2 

sin  B I  <  / - 2- - - 5= 

Y  U2/R2>(1  "  VV  +  l1  "  <Vl/J 


(52) 


If  R2  ■  R^,  then  any  nonzero  B  is  suboptimal.  If  R2  >  so  that 
measurement  1  is  more  accurate,  then  with  A^/A^  “  1  Ecl*  (52)  indicates 
maasuremant  1  is  optimal  for  any  g.  For  a  nontrivial  example,  let 
Rj  ■  2Rjl,  A^  ■  3A2  and  A2/R2  “  1  •  rh,n  E<3*  (52)  indicates  that  6 
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can  b«  as  such  aa  34,5°  before  Che  effectiveness  of  measurement  1 
bacoaas  lass  chan  msasuraaant  2.  Notice  that  for  these  sane 
paramstar  values,  Xq.  (51)  overestimates  the  6  limit  aa  45°.  Evan 
so,  the  noise-level/cosine-squared  relationship  la  a  fair  rule-of- 
thumb  for  eccentric  error  ellipses.  It  should  predict  the  correct 
trends  in  higher  dimensional  error  hyperellipsoida  as  well. 

When  the  eccentricities  are  small  (nearly  equal  uncertainty  in  all 
directions,  then  the  best  choice  is  the  most  accurate  measurement. 

4.3.2  Vector  Measurements 

The  vector  measurement  case  of  Eq.  (1)  is  now  considered. 

In  this  esse  the  noise  coverience  R  is  e  matrix  rather  than  a 
scalar.  Hovevar,  if  R  Is  diagonal  it  is  known  that  the  a  posteriori 
covariance  matrix  of  Eq.  (13)  can  be  computed  by  cycling  through 
Eq.  (38)  once  for  each -measurement ,  The  rows  h  of  H  are  used, 
one  at  a  time  along  with  tha  i1-  diagonal  term  of  R.  When  R  is 
not  diagonal,  an  orthogonal  transformation  on  Eq.  (1)  can  be  used 
to  give 


Az  -  AHx  +  An 

(33) 

or 

z'  “  H'x  +  n' 

(54) 

Tha  modified  measurement  noise  n'  has  a  covariance  matrix 

E{ n'n,T)  -  E{ An  nT  A1} 

■  ARAT  -  R'  (55) 

It  is  always  possiblt  to  sslecc  the  matrix  A  so  that  R’  is  diagonal. 
Than  Eq.  (54)  applies  in  plsce  of  Eq.  (1),  snd  the  noise  covariance 
is  sgsln  diagonal.  For  this  reason  the  assumption  of  a  diagonal  R 
matrix  is  mads  throughout  the  rest  of  the  eigenvalue-eigenvector 
analysis . 

Let  the  eigenvalues  of  Pq  be  numbered  such  that 
^'1  *  X2  >  X3  >  "*  >  ^n’  ^•*u,ne  that  R^  istisfiss  Eq.  (48)  for  all  J, 
se  that  h^  is  tha  bast  single  measurement.  After  processing  h^,  the 
eigsnveluos  of  the  new  coverience  matrix  will  be 

X1  Rl 

X  +  R’  X2*  V  Xn  (56) 
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Of  this  sac,  *2  will  nov  be  the  largest  If 

*1  *2 

*i <  j-rr  (57) 

A1  A2 

In  fact)  if  all  R^  are  equal  so  that  Eq.  (48)  Is  automatically 

satisfied  whenever  X,  >  X  ,  then  the  optimal  measurement  set  would 
T  1  J 

consist  of  h^  ■  j  provided  that  each  R^  satisfies 

X  X 

K  <  Tlf-  <58> 

1  Ai  An 

Eq.  (58)  would  assure  that  each  eigenvalue  in  turn  is  reduced  to 
a  value  leas  than  the  original  smallest  one,  Xft.  If  these  conditions 


were  all  met,  the  final  would  satisfy 


TR(P1) 


Q  X  R 

y  -1- 

i*l  ^i+  ^i 


Eq.  (59)  gives  the  smallest  possible  TR(P^)  for  the  case  where  all 
are  equal,  and  thus  provides  a  limiting  bound  on  measurement 
performance. 

In  actual  fact,  conditions  (57)  or  (58)  are  often  not  satisfied, 
meaning  that  a  smaller  final  error  could  be  achieved  by  using  a 
given  measurement  direction  more  than  once. 

As  a  simple  example,  let 


and  assume  that  three  measurement  possibilities  exist: 


Hi-11  0] 

S2  -  [o  i] 

hj  -  £-1  0] 


.  R]_  -  2 
,  R2  -  1 
.  R3  -  2 


Both  h^  and  hA  are  in  the  direction  of  eigenvector  £^(X^  ■  10) , 
while  h^  is  in  the  direction  of  £2^2  "  E<5*  (^8)  indicates 

that  h^(or  the  best  single  measurement.  If  it  is  used,  the 

first  eigenvalue  decreases  to  5/3  (See  Expression  (56)).  This  is 


'n-pjywrjrrr-r 


not  less  than  end  Eq.  (57)  is  not  satisfied.  The  best  second 

measurement  is  to  use  again  (or  if  questions  about  independent 
measurement  noise  are  of  concern,  use  h^  which  is  the  opposite 
direction.)  Skipping  some  details,  the  final  matrix  for  two 
sets  of  measurements  are 

0 

1 

GDOP  “  ® 

Noise-weighted  GDOP  “  «* 

TR(P1)  -  13/6  -  2.1667 

GDOP  -  2 

Noise-weighted  GDOP  *  4 

This  example  shows  that  a  13X  smaller  error  is  achieved  by  picking 
the  poorest  GDOP  and  Che  noisiest  measurements.  The  noise-weighted 
GDOP  is  sIbo  infinite  for  the  optimal  measurement  set.  This 
example  was  constructed  to  make  the  point  about  the  inability  of 
GDOP  and  noise-weighted  GDOP  to  select  the  best  measurements.  But 
it  is  not  that  an  unusual  an  example,  and  it  is  easy  to  construct 
many  others  of  similar  behavior.  The  same  phenomenon  occurred  by 
itself  when  numerically  processing  GPS-type  four  dimensional  data. 
That  is,  the  worst  available  GDOP  can  give  the  best  performance. 

It  is  pointed  out  that  in  realistic  cases  with  a  fixed  finite 
number  of  h^  vectors  to  choose  from,  it  is  not  generally  true 
that  the  best  pair  of  measurements,  contains  the  single  best 
measurement  aa  one  of  its  members.  Likewise  the  best  set  of  four 
measurements  need  not  include  the  beat  pair  as  a  subset.  The 
results  in  preceding  paragraphs  are  special  cases  where  all 
measurements  are  in  the  direction  of  the  eigenvectors  of  P  . 
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While  these  ere  optimal  results  they  may  not  be  among  the  available 
options  in  a  real  case. 

One  last  topic  is  illustrated  nicely  using  eigenvalue-eigenvector 
concepts— the  notion  of  diminishing  returns.  In  the  preceding  example 
one  measurement  in  the  direction  reduced  X^  fro*  10  to  5/3,  A 
second  measurement  in  the  seme  direction,  and  with  the  same  R  ■  2 
further  reduced  X^  to  10/11. 

Repeated  application  of  the  same  scalar  measurement  generates 
a  decreasing  saquence  of  eigenvalues  (dimension  of  the  error  ellipse) 
at  the  kC^  iteration 


(k*l) 

i 


(60) 


These  satisfy 


0  <  X 


(k+l) 

i 


X 


(o) 

i 


If  these  inequalities  are  normalized  by  the  initial 
obtains  (k+l) 

°-~x7no  - 1 

(k+l)  (k) 

A  parfect  measurement  gives  X  /X  ■  0  while  a 
gives  this  ratio  a  value  of  one.  If  initially  3 


(61) 

,  ,<o) 

value  X^^  ,  one 

(62) 


useless  one 
X<o)  ,  then 


X^/X^ 


x(2)/x  (1) 

a  i  ;  a  1 


2 

3  ’ 


x(3)/xj2>  -2 

i  i  4 


,  (k).  . (k-1) 
Ai  /Ai 


k 

k+l 


Thia  sequence  rapidly  approachaa  unity,  indicating  that  further 
measurements  in  that  direction  are  of  little  or  no  value.  The 
error  ellipsoid  has  been  flattened  to  almost  zero  in  that  coordinate 
direction.  Actually,  thia  daaonstration  of  diminishing  returns 
applies  esactly  to  scalar  maaauraents  in  the  direction  of  an  eigenvector. 
However,  a  similar  effect  is  actually  observed  in  higher  dimensions 
ss  wall.  All  other  things  being  equal,  some  sort  of  alternation  of 
measurements  seems  better  than  a  repetition  of  the  same  sat. 
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Using  the  error  ellipsoid  concept  for  heuristic  justification,  the 
beat  strategy  a ease  to  ha  to  squeeze  down  on  the  largest  dimension 
of  the  ellipsoid,  then  another  large  dimension,  and  so  on,  always 
trying  to  deform  the  ellipsoid  into  the  smallest  possible  sphere. 
Tiae-squeoces  of  measurements  are  discussed  further  in  Section  5. 

4.4  Algorithmic  Approach  to  Measurement  Selection 

The  basis  for  this  section  iB  Eq.  (34).  Consider  the  typical 

jGPS  application  where  measurements  to  four  satellites  constitute 

one  measurement  set.  Assume  that  there  are  N  satellites  visible 

from  which  the  best  four  are  to  be  selected.  (N  would  typically 

be  7  or  less).  Define  the  N  x  4  matrix  H  whose  rows  are  the 

c  . 

vectors  associated  with  a  measurement  to  each  satellite.  If  the  i 
satellite  has  azimuth  angle  A  and  elevation  angle  E,  the  row  is 

•  [coaE  cosA  cosE  sinA  sinE  1]  (63) 

Also  define  the  N  x  N  matrix  R  as  the  covariance  matrix  for  all  N 

c 

noise  components  jointly.  Hc  and  Rc  would  have  to  be  given  data, 
along  with  and  the  weighting  matrix  W.  Two  more  symmetric  N  x  N 
matrices  are  computed  from  the  given  data.  They  are 


B  -  H  P  HT  +  R 

CO  c  c 


C  -  H  P  W  P  HT 


(64) 


CO  o  c 

The  score  associated  with  any  number  of  measurements,  from  1  through 
N,  can  be  obtained  by  computations  on  Iq.  (64).  For  example,  if 
all  N  measurements  were  to  be  used  (and  this  would  give  the  smallest 
possible  a  posteriori  error)  then  Eq.  (34)  gives 


Score  -  TR(CB-1)  (65) 

but  then  there  are  no  options  to  rank.  In  general,  fewer  than  N 
measurements  will  be  processed.  To  indicate  the  number  of  measure¬ 
ments,  as  well  as  which  ones,  a  subscripted  S  will  be  used  for  the 
score.  For  example  S,  means  one  measurement  to  satellite  3  is  made, 
S,^  ass ns  a  pair  of  measures,  to  satellites  i  and  j  are  used,  and 
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S^kl  auni  four  auiurnanti  are  used,  to  satellites  i,j,k  and  1 . 
These  various  scores  are  computed  as  in  Eq.  (65),  but  subsets  of 


tha  matrices  C  and  B  are  used. 


4.4.1  Single  Measurement  Scores 

To  determine  the  relative  ranking  of  N  possible  single  measure¬ 
ments,  the  following  N  scores  must  be  computed  and  ranked. 


Si  “  Cii/Bii  1  "  1,N  (66) 

where  and  are  the  i,i  components  of  C  and  B  of  Eq.  (64). 

4.4.2  Pairs  of  Measurements 

In  the  case  of  measurements  to  satellite  i  and  satellite  j, 
the  2x2  submatrices  formed  from  the  intersections  of  rows  i  and  j 
With  columns  i  and  j  of  B  and  C  are  formed  and  used  in  Eq.  (65) 
instaad  of  the  full  N  x  N  matrices.  The  explicit  form  of  the  2x2 
inverse  is  easily  computed,  giving  the  score 


ij 


^BiiCjJ  +  BJJCii~  2BiJCiJ^BiiBjj“  Bij^ 

n: 


(67) 


There  are  2 ^”2^  “  ^ sets  of  i,J  pairs  that  would  need  to  be 
computed  and  ranked.  The  optimum  i,J  pair  corresponds  to  the  largest 
value  of  It  is  convenient  to  store  the  results  of  Eq.  (66)  along 

the  diagonal  of  an  N  x  N  triangular  matrix  F,  and  the  results  of 
Eq.  (67)  in  the  i , J cBl  off  diagonal  location.  Note  from  Eq.  (67)  that 
S 


a 


Sj^.  Thie  triangular  matrix  F  will  be  used  in  an  approximate 


ranking  scheme  for  four  measurements,  in  Section  4.2.4. 


4.4.3  The  General  Case,  Exact  Formulation 

Although  these  methods  apply  for  any  number  of  measurements 
from  1  to  N,  GPS  is  most  interested  in  4  satellite  fixes.  Attention 
le  therefore  focused  on  evaluation  of  For  each  set  of  ijkl 

values  s  pair  of  4  x  4  matrices  i  and  t  are  formed  from  B  and  C  as 
the  elements  at  the  intersections  of  rows  and  columns  i,J,k,l.  No 
simple  explicit  formulas  are  available  for  the  4x4  inverses  that 
are  required  by 
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(68) 


sijki  "  TR(^  ^"1)  “  ^ 

It  Is  inefficient  to  actually  coapute  the  A  x  4  inverse  and  then 
multiply  by  &  and  then  take  the  trace.  It  is  more  efficient  to  define 
an  unknown  matrix  D  by  the  set  of  simultaneous  equations 

$  D  -  £  (69) 

and  then  use  a  linear  equation  solving  routine  such  as  gaussian 
elimination  to  solve  for  D.  Actually  only  the  diagonal  elements  of 
D  need  be  calculated  since 

Sljkl  ‘  'nl(D>  '  l  °aa  (70> 

J  a*l 

This  is  the  approach  used  in  the  numerical  studies  of  Section  4.2.5. 
An  alternate  method  of  computation  is  available,  using  BAcher's 
formula^  . 

TR(£  If1)  m  -  ~1  TR(2)  -  —  TR(fr  l) 
a4  a4 

-  —  tr(£  l2}-  —  it3)  (7i) 

a4  a4 

where 

ax  •  -TR(fc) 

a2  -  -  y^TR^)  +  TR(&2)J 

®3  "  “  t[*2TR(^  +  aiTR(fi2>  +  ra(^3)] 

*4  *  "  J[«3TR(2()  +  a2TR(Jl2)  +  a1TR(^3)  +  TR(^)]  (72) 

While  this  algorithm  appears  simple,  the  necessity  of  computing 
the  various  matrix  powers  and  products  may  actually  make  it  less 
efficient  than  Eq.  (70).  No  detailed  comparison  of  these  methods 
was  carried  out  during  this  study. 

4,4.4  Approximation  Methods 

Because  of  the  computstional  burden  involved  in  evaluating 
the  exact  score  for  a  large  number  of  candidate  measurement 
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combinations ,  the  desire  to  find  simple  approximations  is  quite 
strong.  Traditional  GDOP  and  noiae-vsightsd  GDOP  can  rightfully 
be  cenaidered  as  examples  of  approximations  to  vhat  Is  really 
wanted.  Various  other  approximations  have  been  investigated 
during  this  study.  One  series  of  approximations  grows  out  of  the 
ine^ua  lities^^ 

TR(£iA  j  (If1)  <  TR(fc  r1)  <  TR(&)1  (r1)  (73) 


min 


max 


or 


r1) 


X  $) 

max 


'.in*4’ 


(74) 


These  inequalities  provide  bounds  on  the  true  score,  and  depend 
on  finding  the  maximum  and  minimum  eigenvalues  of  ^  rather  than 
requiring  that  be  computed.  The  bounds  provided  by  Eq.  (73) 
and  (74)  are  often  so  loose  as  to  be  of  little  value  in  properly 
ranking  measurements.  In  addition,  the  determination  of  eigenvalues 
is  in  itself  a  computational  nuiaance.  In  an  attempt  to  eliminate 
the  eigenvalue  dependency,  another  inequality  was  considered, 


n  X  <  TR(ij)  <  n  X 

min  -  -  max 


(75) 


Eq.  (.75)  combined  with  Eq.  (74)  suggests  chat  perhaps 
TR(2!)/TR($) 


might  have  merit  as  an  extremely  simple  to  compute  method  of  ranking 
measurements.  Numerical  experimentation  along  these  lines  generally 
proved  disappointing. 

It  was  known  from  the  outset  that  the  optimum  set  of  four 

measurements,  ssy  to  satellites  1, 3 , k,  and  1,  need  not  contain  the 

best  pair  of  measurements,  which  might  involve  satellites  i  and  m 

for  example.  However,  it  was  felt  that  all  six  pairs  that  can  be 

mads  from  four  indices  i,J,k,l  ought  to  give  pretty  good  pair  scores. 

That  is,  if  S  is  the  largest  of  all  possible  4-measurement 

scores,  then  S..,  S,.  ,  S,,(  S.,  ,  S,,  and  S.  ,  should  have  reasonably 
ij  Ik*  il*  Jk’  jl  kl 

high  scores  when  compared  with  all  pair  scores.  Because  of  the  ease 
with  which  all  and  scores  can  be  exactly  computed  and  stored 
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in  the  triangular  matrix  F  of  Section  4.2.2,  it  would  be  attractive 

to  approximate  in  tanas  of  elements  of  P.  Considerable  time 

warn  a pent  in  trying  to  discover  a  good  approximation  of  thiB  type. 

Theoretical,  empirical  and  ad  hoc  approaches  to  the  problem  were 

f  9] 

triad.  The  beat  result,  whose  derivation  can  be  found  in1  J  ia 


s2 

ijkl 


■  4  *  4  *  4  +  4  *  4 

+  4  -  -  s2  ♦  s2  ♦  s2] 


where  S  ^  indicates  an  "approximation  to"  This  approximation 

will  be  tested,  along  with  the  exact  result  and  GDOP  results,  in  the 


next  Section. 


4.5  Numerical  Examples  of  Measurement  Selection 


This  Section  evaluates  GDOP,  noise-weighted  GDOP,  the  exact 
score  (with  or  without  the  weighting  matrix  W)  and  the  approximate 

A 

score  S.  The  evaluation  consists  of  seeing  how  well  each  of  these 
methods  is  able  to  rank  35  potential  measurement  combinations  of  4  out 
of  seven  satellites  with  realistic  GPS  geometries. 

The  first  task  is  selection  of  input  data.  The  geometry  (H 
matrices)  were  abundantly  available.  However,  the  measurement 

noise  covariance  matrices  R  were  not.  The  more  usual  way  of 
specifying  measurement  noise  is  in  terms  of  jamming  power  to  signal 
power  (J/S)  ratios,  in  decibel  units. 

The  problem  of  converting  J/S  values  to  noise  covariance 
values  for  the  R  matrix  is  nontrivial  and  depends  upon  the  type  of 
jamming,  the  receiver  characteristics  (such  as  bandwidth,  etc,).  By 
assuming  vide-band  Jamming  noise,  (essentially  white  noise)  and 
narrow-band  signals,  along  with  other  drastic  over-simplifications, 
it  appears  that  the  variance  in  a  given  signal  is  related  to  J/S 
according  to 

Rlt  -  n[l  +  B  10<J/S)*  10]  (77) 

whare  a  is  the  nonjansning  variance  and  8  is  a  function  of  the  receiver 
bandwidth.  Many  J/S  values  seem  to  be  between  about  25  and  90  db. 
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Values  of  a  ■  2.5  and  &  ■  .01  were  used,  as  well  as  a  ■  25,  6  *  .005. 

T  -1  -1 

The  seal*  factor  a  is  not  needed,  if  only  ranking 6  of  (H  R  H)  are 
being  made.  However,  the  more  general  Eq .  (13)  or  (34)  require  the 
scale  factor  be  known  relative  to  values  from  P  .  The  off-diagonal 
entries  of  the  R  matrix  are  related  to  the  angular  separation 
between  measurement  directions,  and  the  antenna  beam  width.  The 
correlation  between  two  measurements  was  programmed  as 


1  2 

Pij  -  2  cos  6ij 


if  e  <  4o 


■  0  otherwise  (78) 

A  few  cases  of  nonpositive  definite  R  matrices  have  cropped  up  using 
this  model.  It  is  not  implied  that  this  noise  model  is  optimum  or 
even  totally  correct.  It  is  felt  that  exact  numbers  are  not 
crucial  to  the  main  points  of  this  study.  When  better  noise  models 
become  available,  they  can  be  used  with  the  selection  methods 
developed  here. 

The  other  key  input  data  is  Pq,  a  4  x  4  matrix.  In  order  to 

proceed  with  the  study  a  limited  parametric  approach  was  taken 

regarding  the  size  of  the  eigenvalues  1  of  Pq.  Values  were  grouped 

into  the  classifications  large,  typical  and  small.  The  eccentricity 

of  the  P  ellipsoid  was  varied  also,  in  terms  of  ratios  of  X./X.. 

o  i  j 

The  values  ranged  from  highly  eccentric  to  spherical.  A  complete 
parametric  study  was  not  carried  out.  Six  different  Po  matrices 
were  tested,  along  with  8  different  sets  of  satellite  geometries, 
several  sets  of  J/S  values  in  some  cases,  and  a  few  tests  with  a  non¬ 
unit  W  matrix.  About  one  hundred  combinations  were  tested.  Eleven 

complete  sets  of  results,  each  consisting  of  two  pages  cf  computer 

[121 

printout,  were  presented  and  discussed  in  .  Some  of  these  results 

[91 

are  also  discussed  in1  ‘ .  Table  1  summarizes  the  key  features  of 
the  numerical  experiments.  Case  numbers  shown  agree  with  earlier 
reports.  In  Table  1  rT  indicates  the  true  rank,  based  on  TR(WP, )  of 
the  measurement  combination  chosen  as  best  by  the  various  methods. 

The  value  rQ  is  the  rank,  as  determined  by  the  various  methods,  of 


33 


Casa 


Brief 

Leecti.p, 


Typical  Pq 
variable  J/S 
W-I 

Large,  more 

eccentric  P 

'  '  —  o 

variable  J/S 
V-I 

Small  P 

o 

variable  J/S 
W-I 

Large  Pq 
equal  J/S 
W-I 

Same  as  case 
5  except 
increased 
correlations 
in  V 

_ o _ 

Typical 
sized  P  , 

C 

but  spheric, 
variable  J/S 
W-I _ 

Small,  highly 
eccentric  P 

o 

variable  J/S 

W-I 


CLOP 

H.W. 

-  CHOP 

1  S 

S 

or 

-TR<P1) 

1  s 

or 

TJKWP^) 

29,16 

2,11 

!■  • 

3,5 

1,1 

1,1 

545.8  . 

412.6 

i  413.9 

410.4 

410.4 

33 

.6 

,  .95 

0  . . . 

0 

3,33 

6,33 

4,2 

1,1 

1,1 

1719 

1747 

1720 

1704 

1704 

.9 

2.5 

.9 

0 

0 

5,35 

2,13 

1,1 

1,1  i 

J- ,  -*• 

41.4 

40.5 

39.2 

39.2  | 

39.2 

5.6 

3.3 

0 

I 

0 

0 

1,1 

2,2 

13,21 

1,1 

1,1 

361.8 

362.9 

677.5 

361.8 

■ 

361.8 

0 

.3 

87 

0 

0 

28,15 

2,8 

3,6 

1,1 

1,1 

455 

280 

282 

279 

279 

63 

.5 

1 

0 

0 

16,1 1 
452 
20 


22,17 

281 

43,7 


1,1 

195.7 

0 


1,1 

195.7 

0 


r„,,r 
T’  o 


Vro 


r_,r 
T  o 


r_,  r 
T  o 


rT’ro 


r_,r 

T  o 


r_  ,r 
T  o 


Tabla  1)  Summary  of  STATIC  Measurement  Selection  lixperiments  (continued) 


Case 


£ 

« 

Brief 

Descrip. 

GDOP 

N.'V. 

GDOP 

A 

s 

S 

or 

TR^) 

s . 

or 

tr(wpl) 

60 

Sum  at 

19JQH 

on 

2,7 

1,1 

1,1 

r„,,r 

T  o 

csss  5  but 

■ 

707 

698 

698 

t 

- 

Wtt-2.5,  1., 
.25,  .25 

46 

m 

1.3 

0 

0 

X 

61 

Same  as 

21,9 

i,i 

6,3 

4,2 

1,1 

Vro 

case  5  buc 

357 

320.4 

327.7 

326.6 

320.4 

e 

W^.25,  2.5, 
.25,  1 

11.6 

0 

2.3 

1.9 

0 

X 

62 

Same  as 

19,16 

2,11 

3,5 

1,1 

1,1 

T  o 

case  5  buc 

309 

259 

259.2 

258.9 

258.9 

e 

1 

WijL-.25,  .25, 
2.5,  1 

19.3 

0 

.1 

0 

0 

*4 

63 

Same  as 

33,16 

4,11 

2,5 

1,1 

1,1 

r^,r 

T’  o 

case  5  but 

461 

333.8 

332.8 

329.3 

325.3 

e 

_ 4 

W11-l,  .25, 
.25,  2.5 

40 

1.4 

1 

0 

0 

% 

the  true  optimum  choice.  For  example,  in  Caae  5  GDOP  selected  as 
baat  a  combination  which  was  actually  only  291*1  best  out  of  35. 

GDOP  placed  the  true  optimum  in  16^  place.  The  value*  of  and  r 
are  all  one  in  the  TR(WP^)  column,  because  the  top  choice  of 
minimum  TR(W?, )  or  majciaua  S  is  by  definition  the  over-all  optimum. 

The  values  e  in  Table  1  ere  the  results  of  TR(WP^)  when  the  top 

choice  measurement  combination  is  used.  In  Case  5,  GDOP  selected  a 

measurement  set  which  was  really  only  29C  beat,  and  this  gave  a 

value  of  TR(WP^)  -  545.8,  compared  with  the  optimal  choice  which 

gave  410.4.  Thus  using  GDOP  caused  a  33%  larger  TRCWP^)  than  the  optimal. 

All  percentages  above  optimal  are  given  in  Table  1, 

Note  that  Case  20  is  an  example  where  the  true  optimal  choice  has 
the  worst  of  the  35  GDOP  values  for  this  geometry.  Since  the  a  priori 

arrors  are  very  small,  as  evidence  by  the  small  P  ,  there  is  little 

o 

to  be  gained  from  the  measurements.  The  absolute  differences  of  the 
reeulta  ara  therefore  small  no  matter  which  selection  method  is  used. 

Case  24  has  a  large  ?g  and  all  J/S  values  are  equal.  This  is 
the  situation  where  GDOP  is  expected  to  do  well,  and  it  does.  The 
noise-weighted  GDOP  does  slightly  worse  even  though  all  diagonal 
terms  in  R  are  equal.  This  is  due  to  measurement  correlations, 
which  cause  nonzero  off  diagonal  terms.  Using  R  ^  as  a  weighting 
term  hurt  more  than  it  helped  in  this  case,  but  the  differences  are 

A 

small.  Casa  24  also  illustrates  the  occasional  example  where  S  is 
not  a  good  discriminator  of  measurement  effectiveness.  Cases  60 
through  64  use  the  seme  input  data  as  Case  5,  except  various  W 
weights  ara  used.  The  trace  of  W  Is  held  at  4  in  all  cases. 

Some  generalizations  can  be  stated  about  the  above  results. 

Tha  use  of  GDOP  to  select  measurement  combinations  involves  com¬ 
puting  the  inverse  of  a  4  x  4  matrix  (or  the  equivalent  computation) 
for  each  potential  measurement  set.  This  procedure  generally  does 

not  give  the  best  selection.  It  becomes  better  as  P  becomes  large 

o 

and  when  the  measurement  noise  levels  are  nearly  equal. 

Adding  the  inverse  of  the  measurement  noise  covariance  to 
GDOP  gives  whet  ie  called  here  noiee-velghted  GDOP.  It  usually 
givea  batter  reeulta  than  GDOP,  if  P  is  large  and  if  the  measurement 
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ooiee  lavala  are  unequal.  There  are  examples  where  this  generaliza¬ 
tion  gives  a  vorae  selection  than  G1>0P  (Ceaea  11  and  20), 

k  reasonably  staple  approximate  ranking  scheme  uses  exact 
analysis  of  measurement  pairs  to  predict  effectiveness  of  seta 

A  ,  - 

of  four  measurements.  This  is  the  S  approximation  method.  It 
often  works  quits  wall,  almost  always  batter  than  CDOP.  The 
occasional  failures  such  as  Case  24,  were  fairly  dramatic,  and  are 
sufficient  reason  to  consider  it  unreliable. 

The  exact  analysis  of  the  weighted  trace  of  the  a  posteriori 
state  covariance  is  the  true  optimal  method.  Its  computational 
burden  is  not  really  much  higher  than  the  4x4  matrix  Inversion 
required  by  the  GDOP  computation.  Several  methods  of  carrying  out 
this  exact  computation  were  discussed  in  Section  4.2.3. 

The  preceding  analysis  and  results  concludes  the  consideration 
of  the  static  measurement  asltction  problem.  It  is  possible  to 
construct  scenarios  whan  the  over-all  best  time  sequence  of  measure¬ 
ments  is  not  the  one  which  results  from  selecting  the  best  set  at 
each  possible  decision  time.  For  example,  if  the  time  needed  to 
acquire  and  settle  on  a  new  combination  of  satellites  is  sufficiently 
long  some  measurement  cycles  may  be  missed.  The  measurement  omissions 
may  contribute  more  to  the  time  average  error  than  would  use  of 
a  slightly  suboptimal  measurement  set,  without  omissions.  Some  of 
these  issues  are  considered  in  the  next  Section. 


5.  SYSTEM  DYNAMICS  AND  TIME- SEQUENTIAL  MEASUREMENTS 

Several  versions  of  the  optimal  measurement  selection  problem 
have  been  considered  in  the  literature,  14, 15, 16]  Alc^OUg^  this 

previous  work  may  have  indirectly  influenced  this  section,  it  does 
not  appear  that  a  ready  solution  to  the  GPS  measurement  selection 
problem  is  yet  available. 

This  section  assumes  that  a  generalization  to  GDOP  is  available 
which  indicates  at  any  given  time  which  combination  of  measurements 
will  yield  the  smallest  mean  square  position  error.  The  question 
addressed  here  is:  What  is  the  optimal  strategy  for  selecting  a 
sequence  of  measurements  over  some  period  of  time?  In  particular, 
does  the  optimal  sequence  consist  of  selecting  the  Instantaneous 
best.  In  the  generalized  GDOP  sense,  at  each  decision  time?  This 
is  a  very  complicated  question  whose  answer  depends  upon  many 
parameters,  conditions  and  ground  rules  as  yet  unspecified.  To 
begin  the  investigation,  a  very  simple  model  is  first  considered. 

5.1  Second-Order  Model 

Let  the  system  state  vector  x  have  Just  two  components, 
position  and  velocity.  Assume  that  at  each  measurement  time  only 
two  choices  are  available: 


*!_  ■  [1  0]  *  + 
or 

z2  -  [0  1]  *  +  n2 


(79) 


Since  z^  is  a  direct  measurement  of  position,  it  will  be  more  effective 
in  reducing  position  uncerteintlss. 

Likewise  z ought  to  be  more  effective  in  reducing  velocity 
uncertainties.  It  will  bscome  clear  soon  that  both  z^  and  z2  affect 
both  position  and  velocity  somewhat,  unless  the  correlation  between 


position  and  velocity  is  zero. 

Let  the  state  covariance  matrix  before  a  measurement  be 


(80) 


Then  the  covariance  P(k|k)  Just  after  a  measurement  la  computed 
from  Eq.  (13)  or  (24).  Let  the  selection  of  measurement  z^  be 
option  A.  Then,  H  ■  [l  oj  and  R  ■  R^ 


A 

f  (k|  k) 


PV 


JPV 


/ 


(i.+ 


Op/^) 


(81) 


Let  the  choice  of  option  B. 

yielding 


B 

P(k|k) 


a 


PV 


Then  H  -  [o 


>1 


and  R 


4/R2} 


'2* 


(82) 


If  the  "best"  choice  is  defined  as  the  one  which  gives  the  smallest 
position  uncertainty  after  just  one  measurement,  then  a  comparison 


A  B 

of  the  1,1  elements  of  P(k|k)  and  P(k|k)  is  appropriate.  The 
choice  depends  on  several  parameters,  even  in  this  simple  case. 

But  as  p  o  option  A  clearly  is  the  best  choice.  This  direct 
measurement  of  position  would  be  expected  to  be  the  best  way  to 
reduce  position  uncertainties  in  general,  (although  an  unusual 
combination  of  p,  oy/ R2  and  Op/R^  could  cause  this  not  to  be  true). 

Vhen  the  "bast"  choice  is  defined  to  mean  the  one  which  gives 
the  smallest  average  position  error  over  Che  interval  T  between  two 
measurement  times,  the  decision  is  more  interesting.  To  investigate 
this  question,  the  transition  matrix  $  is  needed,  as  will  as  the 
plant  noise  covariance  matrix  Q.  A  reasonable  approximation  for  $ 
in  many  cases  ia 

fl  T  1 


♦ 


1 
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Than  the  extrapolated  covariance  matrix  (after  making  a  measure¬ 
ment  and  than  extrapolating  up  to  just  prior  to  the  next  measurement) 
ia  P(k+l|k)  computed  as  in  Eq.  (26).  If  option  A  ia  selected,  the 
1,1  component  ia 


A 

PCk+l|l0lfl 


{Op  4-  2T  opv+  ojfl  4-  Op(l-p2)/RjT2} 

1  +  4/Rl 


(83) 


If  option  B  is  selected,  the  1,1  component  is 

B  {o2  [l  +  O2(l-02)/R  1  4-  2T  o  +  o2  T2} 

P(k+l|k).  . - - - - - \ - — - - -  (84) 

1  +  c2/R2 

These  two  expressions  could  be  compared  if  the  desire  is  to  get  the 
smallest  extrapolated  position  arror.  However,  it  seems  more 
appropriate  to  seek  the  smallest  average  error  over  the  interval  T. 

If  a  straight  Hue  variation  ia  assumed  between  the  beginning  value 
(from  Eq.  (81)  or  (82))  and  the  ending  value  (from  Eq.  (83)  or  (84)) 
then  the  appropriate  thing  to  do  is  compare  the  averages 


A  A  B  B 

\  l[p(k|k)  4-  P(k+l|  k)J  >  j  T  [p(k|k)  +  P(k+l|k)]  (85) 

Carrying  out  the  algebra  and  subtracting  side  B  from  side  A  gives 
the  decision  variable  D.  Mote  that  the  subtraction  causes  the 
plant  noise  terms  to  cancel. 


°  ’  [2fe5)2  +  2(o)  »  +  l]  <»*'  *i 

V  V 


4-  0;(i-pfc)|^  - 


r2t 


(86) 


where  g^  -  1/(1  4-  Op/R^)  end  gg  -  1/(1  4-  Oy/R^).  D  <  0  means 
choice  A  is  best,  while  D  >  0  means  choice  B  is  best,  and  D  *  0 
means  both  choices  are  equally  good. 
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The  trade-off  can  beat  be  described  in  terms  of  Figure  3. 


Figure  3.  Linear  Error  Growth 

Choice  A  will  normally  give  a  smaller  initial  position  error  but  a 
large  velocity  error.  The  velocity  error  determines,  at  least  in 
part,  the  slope  of  the  error  growth  between  measurements.  Even 
for  the  simple  problem  under  consideration,  there  are  six  parameters 
that  go  into  the  decision,  namely  cp,  o„,  o,  T,  R.  and  R  . 

f  v  2  2 L  1 

From  Eq.  (86)  it  is  seen  that  if  Op/R^  “  Oy/^  (sort  of  a 

signal  to  noise  ratio  SNR,  then  g.  -  g_,  If  in  addition  o2  -  1,  or 
2  Ad 

if  R  -  R,T  / 2,  chen  D  -  0.  Combining  the  two  conditions  not 

1  ^  2  2  2 
involving  p  shows  that  D  -  0  if  the  SNR  are  equal  and  if  o  T  -  2o  . 

2  2  2  v  f 

If  ov  T  >  20p  then  option  B  is  best. 

There  are  many  combinations  of  parameters  that  lead  to  choice  B 

being  better  than  choice  A.  This  suggests  a  partial  answer  to  the 

question  which  opened  this  Section.  The  choice  of  measurements 

which  gives  the  smallest  instsntaneoua  position  error  is  not 

necessarily  the  one  which  will  give  the  smallest  time-average 

[13] 

position  error.  A  very  similar  system  model  has  been  considered1 
and  it  is  shown  that  if  only  one  type  of  measurement  can  be  used, 
the  position  measurements  are  optimal.  When  either  position  or 

velocity  measurement  can  be  used,  it  is  found  that  a  combination 

[13] 

of  both  is  best  over  a  period  of  time.  These  results1  apply  for 
a  specific  set  of  parameters  equivalent  to  the  six  mentioned  in 
the  preceding  paragraph.  An  examination  of  a  more  realistic  model 
of  the  GPS  system  is  now  presented  to  show  that  similar  conclusions 
also  apply  there. 
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5.2  Mora  Gmnartl  GPS  Modal 

Let  the  state  vector  x  consist  of  seven  components;  three 
position  components  plus  a  time  bias  fora  the  first  four  components, 
and  these  four  have  been  dealt  with  extensively  in  earlier  sections. 
The  last  three  elements  of  x  are  the  user  velocity  components. 

The  state  covariance  matrix  F  is  7  x  7  end  will  be  partitioned  as 


5  -  I - 1 


where  is  the  4x4  block  dealt  with  in  previous  sections, 


the  3x3  coverlance  of  the  velocity  terms  and  is  the  4x3 
matrix  of  correlations. 

It  is  assumed  that  measurements  are  processed  in  set9  of  four. 
That  Is,  r  is  a  four-component  vector  related  to  the  state  by 

*  -  U  x  +  n  :88) 

The  measurement  matrix  U  ie  4  x  7  and  can  also  be  partitioned  as 

“  ‘  f“l  i  °4x3  ]  <"> 

where  Is  the  4x4  measurement  matrix  of  previous  sections. 

The  measurement  noise  covariance  matrix  R  is  also  of  dimensions  4x4. 

As  the  result  of  processing  the  k+l^  group  of  four  measurements 
the  covariance  is  updated  by  applying  Eq .  (24),  which  in  partitioned 
form  becomes 

P  (k+l|k+l)!p  (k+l|k+l)  f  P  (k+l|k)|P 2(k+l|k) 

-i - 4-f - *  .  i  -  -  -  t  - - 

PUk+lIk+D'.P  (k+ljk+1)  J  P,(k+l|k):P,(k+l|k)  _ 


P1(k+l|k)H^[H1P1(k+l)k)Hj+  R]"1H1P1(k+l|k)  ! 


P2(k+l|k)Hjfli1P1(k+l  |k)Hj->-  r]  “1H1P1  (k+1  |k)  ! 
P^(k+l|k)H^  Jhi P :  (k-t-l  |  k)H^+  R]“1H1P2(k+l|k) 
P2(k+l|k)H2  ^1PL(k+l|k)H2+  (k+1 j k) 


It  is  stan  that  tht  uppar  laft  hand  4x4  partition  of  Eq .  (90) 
is  exactly  tha  same  fora  as  was  dealt  with  vhan  considering 
generalised  GJXJP  performance  measures. 

In  order  to  analyse  time-sequential  behavior  a  state  transition 
matrix  ♦ (k+1 , k )  and  a  plant  noise  covariance  matrix  Q(k)  are  required. 
Aaauaa  that  the  motion  of  the  user  can  be  approximated  as  rectilinear 


motion  over  the 


ill  time  T  ■  -  t^,  and  that  the  time  bias 


remains  constant  except  for  random  effects.  Then  the  7  x  7  $  matrix 
can  be  approximated  by 


x3  ;°i  x3 T 

-  -  -I-  i-  -  -  -  (9] 

0 (k+1  ,k)  -  0  0  0<1  |0  0  0 

!o ; 

.  °3x3io!  13  _ 

where  1^  and  °3x3  are  the  3x3  unit  and  null  matrices  respectively. 

If  all  plant  ncise  effects  on  position  and  velocity  are 

attributed  to  3  independent  components  of  random  acceleration  with 
2 

equal  variances  o  ,  and  if  the  variance  of  the  random  time  bias 
2  2  a 

drift  is  6  o  T  ,  then  the  7  x  7  Q  matrix  can  be  approximated  as 

d 


t4I  0 
U  T  1  0 
3  4  '0 


0  0  0 


0  0  0  0 


Using  these  approximations,  the  covariance  propagation  equation  of 
Eq.  (26)  can  be  written  as 


P(k+1  k) 


r  h]  t 

Jl+  -  -  “  P2 

to  0  0J 


T+  P,  I,.  0  T+  -  P,  I  ;0  T  •  P,+  -  -  -  P ,T 

2 1  3|  0  1 0  0  oJ  3l  3:  oj  :  2  0  0  0  3 


P2  +  P3fI3iolT 


+  Q(k) 


(93) 


The  trace  of  the  upper  left  4  x  4  partition  of  P  will  be 
investigated;  since- this  is  in  essence  what  the  GDOP  and  generalized 
GDOP  approaches  use.  Define 

e(k+l|k)  -  TR  [p x (k.+l |  k)]  (94) 

and  e(k+l|k+l)  -  TR  jp^<k-»-l|  k+l>]  (95) 

Then  from  Eq.  (92)  and  Eq.  (93)  2 

9  3o  ?  ? 

e(k+l  |  k)  -  e(k|  k)  +  2aT  +  bT  +  ~~  T  +  Bo^  T  (96) 

where 


a  »  TRACE  of  the  first  3x3  sub-block  of  P  (k|k) 

b  -  TR[P3(k|k)]  (97) 

From  Eq.  (90) 


e ( k+1 | k+1 )  -  e (k+l| k)  -  S(k+1)  (98) 

where  S(k+1)  is  the  "score"  associated  with  the  measurement  set  selected 
at  time  k+1,  as  introduced  in  Eq.  (33)  and  (34)  but  with  W-I. 

To  be  consistent  with  Eq.  (85)  of  Section  5.1,  potential 
measurements  (defined  by  choice  of  and  R)  could  be  compared  based 
on  the  average 

i  T  Je (k| k)  +  e(k+l|k)]  . 


It  can  be  seen  that  exactly  the  same  kinds  of  trade-offs  exist  here 
as  in  Section  5.1.  That  is,  some  candidate  measurement  set  may  produce 
the  lowest  e(kjk)  (generalized  GDOP)  at  time  k,  but  may  not  be 
optimal  because  of  a  larger  associated  velocity  uncertainty  as 
represented  by  b  in  Eq.  (96)  and  Eq.  (97). 

Instead  of  Just  looking  at  a  single  time  step,  the  average 
position  error  over  K  time-cycles  (the  planning  horizon)  ought  to 
be  minimized  by  choice  of  measurements.  That  is, 

K 

minimize  J  ■  £  e(k+l|k+l)  (99) 

k-o 
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by  considering  All  admissible  measurements >  subject  to  the  recursive 
equation,  obtained  by  comparing  Eq,  (96)  and  (98). 

e(k+l|  k+1)  -  e(kjk)  +'2«kT  +  (bk+  8u2)T2 

3c2  . 

+  — |  T  -  S(k+1)  (100) 

The  decision  variable  is  the  score  S(k+1).  The  formulation  of 
Eq.  (99)  and  (J. 00)  looks  like  a  standard  type  of  discrete  optimal 
control  problem  that  would  be  amenable  to  solution  perhaps  using 
dynamic  programming.  This  is  not  exactly  true,  because 

(1)  the  score  S(k+1)  is  in  fact  a  function  of  e(k+l|  k)  and 
thus  indirectly  e(k|k) 

(2)  the  parameters  ak  and  b^  are  not  constants  but  are  affected 
by  the  choice  of  measurements 

(3)  the  real  crux  of  the  GPS  problem  is  the  specification  of 
the  set  of  admissible  measurements.  Obviously,  satellites 
must  be  above  some  minimum  elevation  angle  if  they  are  to 
be  admissible  candidates.  More  than  this  however,  a 
variety  of  ground  rules  could  be  Imposed.  For  example, 
once  a  given  satellite  is  selected,  it  may  be  required 

to  stay  with  it  for  some  minimum  time.  Also,  to  allow 
for  acquisition  and  settling  times,  a  new  satellite 
may  not  be  admissible  until  after  some  switch-over  delay. 
Thus,  the  admissible  measurements  change  with  time  in 
some  complex  way. 

It  is  possible  that  Eq.  (100)  might  be  made  a  useful  approximate 
model  by  assuming  some  typical  average  values  for  and  b^ 
(correlations  and  velocity  uncertainties)  and  by  approximating  S(k+1) 
by  a  measurement  effectiveness  indicator  analogous  to  the  original 
GPOP  which  does  not  depend  on  the  covariance,  This  type  of 
formulation  could  then  be  solved,  using  dynamic  programming  for 
example.  This  approach  would  have  to  be  investigated  by  means  of 
extensive  realistic  simulations,  and  is  therefore  outside  the  scope 
of  this  study, 
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The  um  example  of  Section  2 .3  is  reconsidered.  This  could 
be  treated  as  the  2x2  position  sub-block  of  a  4  x  4  system  with 
♦  ■  I  and  Q  ■  0.  Alternately  it  might  be  repeated  measurements 
on  a  static  system.  At  each  measurement  time  a  choice  must  be 
■ads  between  measurement  seta  1  and  2.  both  have  the  sane  noise 
matrix  R,  and  the  initial  covariance  is  the  same  one  used  in 
Section  2.3.  In  that  previous  section  it  was  shown  that  set  1  has 
a  slightly  better  GDOP  but  set  2  gives  a  smaller  mean  square  error 
after  one  measurement  update.  Here,  all  eight  sequences  of  measure' 
manta  over  the  first  three  measurement  times  are  analyzed.  Major 
results  ere  given  in  Table  2.  The  initial  mean  square  error  in  all 
eight  cases  is  TR(Pq)  ■  6.  The  successive  values  of  e(l| 1) ,  e (2  j  2) 
and  e(3|3)  defined  in  Eq.  (98)  or  (100)  are  tabulated. 


Measurement 

Sequence 

1st  2nd  3rd 

e  ( 2 1  2 ) 

m 

l  e(i|  i) 
i-1 

H1 

H1 

H, 

X. 

2.823 

1.895 

1.433 

6.151 

*1 

H1 

H2 

2.823 

1 

1.895 

1.352 

6.070 

H1 

H2 

H1 

2.823 

1.907 

1.433 

i 

6.163 

H1 

H2 

U2 

2.823 

1.907 

1.443 

6.173 

H2 

H1 

Hi 

2.697 

1.782 

1.352 

5.331 

H2 

E1 

H2 

2.697 

1.782 

1.352 

| 

5.831 

H2 

H2 

H1 

2.697 

1.856 

1.352 

l 

5,905 

H2 

H2 

H2 

2.697 

1.856 

1.433 

5,936 

Table  2:  A  Sequence  of  Three  Measurement  Updates 


There  are  essentially  just  two  different  "final  values"  for  e(3|3). 
Tha  larger  value  of  1.4+  occurs  when  or  5^  is  used  three  times 
In  a  row.  Tha  smaller  and  value  of  1,35  is  obtained  when  some 
variety  exists  in  choice  of  or  It  is  also  noted  that  the 

sum  (or  average)  of  the  three  e(. |i)  terms  has  a  more  selective 
minimum,  and  it  occurs  for  the  alternating  sequences  H  , 

(although  H^,  H^,  gives  virtually  the  same  result).  This 
demonstrates  that  the  advantages  of  alternating  measurements 
discussed  in  Section  4.3.2  also  apply  in  this  case. 

The  discussion  of  diminishing  returns,  in  conjunction  with 
Eq.  (62)  applies  approximately  to  the  mean  square  error  terms 
here.  That  is, 


6  2’  a a|l) 


2  e(3 1 3)  2 

3  '  e (2  j  2)  3  4 


which  fits  the  pattern  given  earlier  for  the  single  axis 

eigenvalues.  This  hardly  constitutes  a  proof,  but  it  does  suggest 

that  in  setting  up  an  optimal  sequential  selection  problem  (like 

Eq.  (99)  and  Eq.  (100)),  some  sort  of  discounting  factor  should  be 

applied  to  each  successive  use  of  a  measurement  set  with  a  basic 

kSQ 

score  of  Sq,  Perhaps  S(k+1)  ■  would  be  a  reasonable  approxi¬ 

mation.  The  determination  of  Sq  would  depend  upon  the  geometry, 
noise  levels  and  state  uncertainty  existing  at  the  beginning  of 
the  sequence.  The  simple  pattern  demonstrated  here  will  be 
disrupted  somewhat  where  ♦  i  I  and  Q  f  0.  Additional  study  would 
be  needed  to  clarify  this  point. 
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6.  SENSOR  SUBSTITUTION  STRATEGIES 

It  la  vail  known  chat  the  smallest  sun  square  estimation.  ' 
error  results  when  every  available  measurement  is  used.  Even  a 
very  noisy  measurement  helps  a  little,  or  at  least  does  not 
degrade  accuracy,  assuming  that  the  appropriate  noise  covariance 
Is  used.  The  same  statement  applies  to  any  secondary  sensors  such 
an  radar  altimeters.  If  this  philosophy  was  adopted  there  would 
be  no  measurement  selection  problem.  The  procedure  would  be  to 
Just  use  everything. 

Hardware  and  software  limitations  normally  prevent  the 
implementation  of  the  "use  everything"  approach.  For  example  the 
number  of  viaible  satellites  will  often  exceed  the  number  of 
receiver  channel)  available  for  tracking.  Therefore  choices  must 

be  made. 


6.1  An  Optimisation  Problem 

One  approach  to  theae  types  of  problems^’**^’^  iB  :o  define 
an  over-all  maaauremant  vector  consisting  of  very  possible 
measurement.  Each  vector  subset  is  multiplied  by  a  scalar 


IB  Z 

l  -l 


m,  z„ 

.  2  -L 


(101) 


Then,  a  scalar  cost  function  which  depends  upon  the  resulting 
mean  aquare  error  (auch  aa  Eq.  (99)  for  example)  and  the  cost  of 
taking  and  processing  each  measurement  is  minimized  with  respect 
to  a  ,  A  common  approach  is  to  use  the  matrix  version  of  the 

1  f  18] 

minimum  principle1  ,  with  the  elements  of  the  covatiance  matrix  P 

aa  state  variables.  If  a  given  scalar  ma  ■  0,  that  means  the  lCtl 
group  of  measurements  Is  not  used,  while  m^  ■  1  means  it  is  used. 
Different  restrictions  on  m^  have  been  considered,  including 
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I®.  ■  l^^and  a .  ■  0  or  appears  that  the  essence 

of  tha  results  in1  J  are  to  not  include  any  measurement  whose 
coat  ezceeda  lta  value,  and  the  "value"  leeks  vary  similar  to  what 
wee  defined  as  the  score  in  Eq.  (34). 

If  the  optimal  selection  measurement  problem  is  formulated 
as  in  Eq.  (101),  and  if  all  measurements  are  included,  whether 
they  be  from  GPS  satellites  or  from  miscellaneous  alternative 
devices,  then  the  sensor  substitution  problem  would  not  need  to  be 
considered  as  a  separate  iosue.  Rather,  it  could  be  included  with 
topics  of  Sections  4  and  S. 

b. 2  A  Suboptlmal  Approach 

The  approach  taken  here  ie  a  auboptimal  one,  based  on  maintain¬ 
ing  separate  identities  of  GPS  satellite  measurements  and  those 
derived  from  nominally  less  accurate  backup  sensors. 

Let  the  full  state  vector  x  which  ia  being  estimated  be 
partitioned  into  two  parts,  and  let  the.  satellite  measurements  z  also 
be  partitioned  into  two  parts. 


Theu  the  measurements  of  Eq.  (1)  can  be  written  in  partitioned 


_li  .  -iii  1 III  +  li_ 

_  lz  .^21  !  H22.  -2  -2 


(102) 


The  question  under  study  assumes  that  another  completely  separate 
method  of  solving  for  the  part  of  the  state  is  available  from 
alternate  (non-GPS  satellite)  sensors.  That  is,  measurements 
£3  satisfy 


-3  “  H32  -2  +  -3 


with  H32  invertible  so  that 


x  “  H_1  z 
-2  n32  -3 


(103) 


(104) 
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Dm  alteraate 
good  CM  aunriNMeti,  ^ 
CPS  Mnw»j\|a.  Art 
contained  In  the  martlet 
placed  on  the  else  of  s^. 


inferior  to 
thee  degraded  or  heavily  jammed 
aoeecte  which  are  degraded  ere 
o  particular  restrictions  are 
I'M  repreeent  pseudo-range  data 
‘."fieviea  could  he  e  single  ttate 


from  one  or  gore  satellite! 
coapooent  (euoh  ee  eltitad*  trow  a  radar  eltiaeter)  or  aore  then  one 
coaponent  (eueh  ae  the  borlaontal  position  from  an  inertial 
navigator).  It  la  known  that  the  true  optimal  estimate  of  x  would 
use  all  three  sources  of  data,  *&4  *3  (assuming  that  the  noise 

atatiecica  for  the  degraded  vara  known).  The  state  update 
equations  from  Eq.  (12b)  can  be  expressed  in  the  partitioned 
form  as 


A  +  a* 

2l  “  2i 


~  H11  il  “  H12  - 2 ^ 


and 


+  K12(z2  -  H2l  xx  -  H22  x2) 
+  Kl3(-3  "  H32  -l'* 


-2  "  -2  +  K21^1  “  H11  -l  ~  H12  -2^ 
+  K22(z2  -  H21  x1  -  H22  x2> 


(105) 


K23(-3 


H32  x2) 


(106) 


where  superscript  -  and  +  indicated  before  and  after  the  measurement 
update  raapectively .  One  alternate  sensor  strategy  consists  of 
dropping  z2  entirely  and  depending  on  z^  to  estimate  x2>  That  is, 
a  possible  algorithm  form  might  be 


5i 


*+  -I 
-2  "  H32  -3 


-L  *  Kll(5l  "  H11  -l'*  +  K13  -3 


(107) 

(108) 


Equation  (108)  reduces  to  the  same  form  as  Eq.  (105),  if  K 
and  if  K 


13 


K11  H12  H32 ' 


12 


Also,  Eq.  (106)  ia  the  same  as  Eq.  (107) 

11  ^23  ana  approacn  zero  and  K2j  ■  H^.  Questions  of  interest 

are:  Whan  can  z 2  be  dropped  and  still  maintain  near  optimal  solutions? 
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J  . 


Of  several  potential  alternatives  (different  partitions  of  *2), 
which  one  is  the  beat  choice?  The  aethod  of  approaching  these 
Issues  is  to  first  determine  detailed  expressions  for  the 
partitioned  euboptiaal  gain  matrix  Inherent  in  Eq.  (107)  and  (108), 


V 


'*ll 

0 

0 

0 

» 

■  -  - 

itioned 

form 

'*u 

*12 

,K21 

*22 

-*UH12H32 

h"1 

H32 


(109) 


Then  the  partitioned  form  of  the  optimal  gain  matrix 

*13 

*23  J 


(110) 


Is  determined.  This  is  done  by  applying  the  usual  Kalman  gain 
equation 

K  -  P  HT[rPHT  +  R]"1  (111) 

to  the  full  set  of  measurements  consisting  of  Eq.  (102)  plus  (303). 
Then,  conditions  under  which  K  approaches  K1  when  noise  levels  n2  in¬ 
crease,  will  be  presented. 

6 .3  The  Suboptlmal  Gain  Matrix 

The  form  of  the  suboptlmal  gain  matrix  K*  is  completely 
specified  once  Is  determined,  Motivated  by  Eq.  (103)  and  the 
first  part  of  Eq.  (102),  define  an  equivalent  measurement 

2  "  2i  *  Hi2  H32  -3 

-  Hu  xx  +  (n1  -  H12  h"2  n3)  (112) 

If  n^  and  n^  are  independent  and  have  covariance  matrices  R^  and 
R33  respectively,  the  equivalent  total  noise  v  ■  -  H^2  n^ 

-1  -T  T 

has  covariance  +  H^2  H32  R^3  H32  Hj^.  The  optimal  estimate  of 
from  Eq.  (112)  Is 
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-  -4i  *  +..cov.<al,g)  [pov  Cj.a)]  -1  [a  -  E{a}] 

"  *1  +  Kutil  “  Hu  ill  “  *11  H12  H32  -3 
where  ■  covfx^.m]  [covCa.a)]  1 

Using  the  covariance  definitions 

cov^a.b]  ■  E{[a  -  E(a)]  [b  -  E(b)]T} 
immediately  ludi  to 

K11  “  P11  P11  U11  +  R11  +  H12  H32  R33 


(113) 

(114) 


11  P11  U11  +  R11  +  H12  H32  R33 


fT  HX  V1 
32  H12J 


(115) 


6.4  Comparison  of  the  Optimal  and  Suboptlmal  Gains 

The  partitioned  fora  of  the  a  priori  covariance  (Just  before 
the  aeasureaent)  is  ^ 


► 

11 

PI2 

■  )T 
12; 

P22 

The  superscript  minus  signs  will  be  omitted  for  convenience  in 
this  section.  The  full  measurement  matrix  is 

\l  R12 ' 

H  ”  H21  H22 

°  h32 

and  the  measurement  noise  covariance  is 


[  0  0 

Substitution  into  Eq.  (Ill)  gives 


T 

T 

i  T 

,  T 

T 

1 

11 

_H11 

+  ?12 

“12 

1  P  H  +  P 

1  rll  21  12 

X  _  _  _  _  _ 

H22 

1  P  H 

•  r12  h32 

»  _  _  _ 

u'1 

T 

T 

T 

!  _T  T 

,T 

1  ,,T 

) 

12 

H11 

+  P22 

H12 

i  V12  H21  +  V 22 

H22 

!  ^22  h32  _ 

(116) 
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where 


*nhl*h  +  hlhA*  H11PUH22  *  H11P12U22  1  H11P12H32 
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If  the  aubaatrix  partitions  of  U  are  called  for  convenience, 

and  if ,  dv 
neglected 


and  if,  due  to  ^  becking  large,  the  following  term,  cen  be 
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then  Eq.  (116)  and  (119)  yield 
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PllHll(Ull”  U13U33U23) 
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reduce*  to  “KnHi2H32  Provlded  th4t  ui3U33  "  Ui2H32’  K2l  is 
z«ro  if  ^32^33^31  “  ®i2  of  thess  conditions  are  set  if 


R33  **  U32P22H32 


(122) 


Eq.  (122),  plus  a  matrix  inversion  Identity,  was  also  used  in 


-1 


reducing  a  long  expression  for  down  to  just  Finally,  the 

lnvarse  matrix  in  is 
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,  -1  T  -1 

11  "  K11  Provlded  that  F22  <<;  H32R33H32’  Howav*r  this  is 


Just  an  alternate  form  of  Eq.  (122).  Therefore  K  -*  K  provided 
R22  becomes  efficiently  large  so  that  terms  (118)  are  negligible, 
and  Eq.  (120)  and  (122)  are  satisfied. 


6.5  Application  of  the  Substitution  Strategy 

If  a  given  measurement  or  group  of  measurements  becomes 
very  degraded,  so  that  increases  to  a  very  large  value,  then 
an  alternate  method  of  solving  for  one  or  more  states  should  be 
considered.  If  more  than  one  substitution  choice  is  available, 
the  one  which  gives  most  nssrly  uncorrtlated  from  x^  ,  should 
bs  aslscted  eo  that  Eq .  (120)  is  most  nearly  true.  The  substitu¬ 
tion  should  be  made  when  becomes  sufficiently  degraded  (due 
to  the  increasing  ^2^  *°  Eq.  (122)  is  reasonably  satisfied. 


55 


That  is,  under  normal  oparation  it  ia  assumed  that  the  secondary 
aanaor  ia  laaa  accurate  than  GPS  navigation  accuracies,  so 
*33  *  *35*2^32*  dagreded-oonditiona  tha  aanaa  of  tha  . 

inequality  will  evantually  reverie. 

Soma  of  the  ideal  presented  above  are  similar  to  those  given 

previously,  for  other  purposes,  .by  Stubberud  and  Wiener  . 

Tha  requirement  that  the  state  vector  be  partitioned  into  the  moat 
nearly  uncorrelated  subvectors  ^  and  x2  is  one  of  the  major  points 
of  aimilarity.  The  best  method  of  verifying  the  appropriateness 
and  effectiveness  of  this  procedure  is  through  realistic  simulation 

tests. 
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7.  ANALYTICAL  TESTS  AND  SUGGESTED  SIMULATION  CASES 

7.1  Analytical  Testa  -  ~ 

A*  stated  -in  Section  3,  one  of  the  casks  of  this  study  was 
to  carry  out  analytical  testa  of  the  proposed  GDOP  generalizations 
and  extensions.  This  task  has  been  accomplished  in  earlier 
sections  concurrently  with  methodology  development.  Specifically, 
in  Section  4  examples  were  presented  which  demonstrated  that  the 
best  GDOP  does  not  guarantee  highest  accuracy.  It  was  shown  that  the 
eigenvalues  and  eigenvectors  of  the  covariance  matrix  play  -a  very 
important  role,  in  addition  to  geometry  and  noise  levels.  Analytical 
tests  were  given  in  low  dimensional  cases,  and  these  were  used  to 
extract  some  general  principles  and  rules  of  thumb  for  higher 
dimensional  cases. 

In  Section  5  analytical  tests  on  a  low  dimensional  model 
demonstrated  that  the  measurement  types  which  are  best  for  the 
instantaneously  minimum  position  error  need  not  give  the  minimum 
average  position  error  over  time.  A  law  of  diminishing  returns 
for  repeated  use  of  the  same  measurement  was  developed  analytically, 
as  was  the  suggestion  that  a  sequence  of  measurements  will  be  more 
effective  if  there  is  some  variety  in  the  sequence. 

7 . 2  Suggested  Simulation  Cases 

A  considerable  amount  of  computational  work  was  carried  out 
in  conjunction  with  this  study.  Some  of  it  applied  to  simple  examples 
but  much  of  it  applied  to  realistic  GPS  geometry.  The  measurement 
noise  model  used  in  these  cases  was  a  hypothetical  one.  The  a  priori 
state  covariance  matrices  were  selected  parametrically.  It  is 
suggested  that  the  methods  of  Section  4,  and  in  particular 
Saction  4,4,  be  tested  in  a  realistic  full  blown  simulation  with 
the  latest  noise  models  end  covariance  matrices.  Even  the  geometry, 
which  was  chosen  realistically  at  the  time  of  this  study,  deserves 
additional  attention  because  of  the  smeller  number  of  satellites 
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in  the  scaled-down  GPS  network.  This  affects  viewing  tines  and 
reduces  Che  number  of  potential  measurement  choices  at  any  given 
time.  One  type  of  case  that  is  expected  to  benefit  most  from 
generalised  GDOP  methods  has  a  highly  eccentric  state  covariance 
error  ellipsoid,  a  vide  range  of  measurement  noise  levels,  end 
GDOP  values  that  ars  all  reasonable  or  roughly  comparable  with 
each  other. 

In  Section  4.5  it  was  shown  that  root  mean  square  position 
error  could  be  improved  by  e  few  feet  up  to  maybe  10  feet  at  a 
given  measurement  time  by  using  s  generalized  GDOP.  Although 
smell,  these  improvements  ere  significant  in  view  of  GPS 
accuracy  goals.  This  level  of  static  improvement  should  be  verified 
with  a  full  and  realistic  simulation.  Such  a  simulation  could  also 
determine  how  these  relatively  small  improvements  in  accuracy  propagate 
in  time.  Do  they  remain  smell  or  do  they  accumulate,  and  if  so, 
how  much? 

The  entire  subject  of  sequential  measurements  needs  to  be 
studied  by  means  of  simulation.  One  scalarized  version  of  the 
problem  was  formulated  in  Section  5.2,  Some  remaining  details  of 
thet  method  can  only  be  extracted  from  actual  simulations.  If 
that  scalar  approach  to  setting  up  the  problem  proves  feasible, 
then  a  discrete  time  dynamic  programming  solution  over  a  finite 
planning  horizon  can  probably  be  obtained  in  near  reel  time.  If 
the  scalar  approach  is  not  satisfactory,  then  the  higher  dimensional 
full"  matrix  approaches  of  References  [13]  through  [16]  could  be 
considered.  The  additional  complexity  would  be  a  disadvantage  for 
real  time  applications. 

A  reasonably  simple  sensor  substitution  strategy  was  detailed 
in  Section  6.  This  should  also  be  evaluated  by  simulating  one  or 
more  very  bad  pseudo-range  measurement  channels.  The  partitioning 
of  the  state  into  essentially  uncorrelated  subparts  is  a  general 
requirement  that  can  only  be  tried  in  simulation.  Also,  several 
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Mtrlx  inequality  condition*  were  developed  in  Section  6.4. 

Exactly,  whet  it  meant  by  one  matrix  being  amall  compared  to  another 
va,  not  stated.  Several  measures  of  relative  emallnesi  can  "be 
considered,  but  only  simulation  can  give  a  feeling  for  which  one 

isorks  best. 


5V 


8.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  study  has  shown  chit  measurements  selected  because  - 
thsy  hsvs  ths  beat  GDOP  do  not  necessarily  give  the  smallest 
estimation  error.  In  fact  examples  have  been  presented  in  which 
the  worst  GEDOP  measurement*  give  the  best  position  estimate. 

Generalisations  to  GDOP  have  been  developed.  Some  of  these 
are  approximations,  and  ona  la  exact.  Performance  improves, 
but  at  the  expense  of  requiring  more  input  date  and  a  higher  computa¬ 
tional  burden.  The  ability  to  minimise  a  weighted  sum  of  the  : 
estimation  error  components  is  built  into  the  exact  algorithm  and 
this  may  be  useful  in  some  applications.  In  addition  to  the 
development  of  an  exact  algorithm  for  ranking  measurement  effective¬ 
ness,  eigenvalue-eigenvector  theory  was  used  to  obtain  insight 
into  the  tradeoffs  among  geometry,  noise  levels  and  a  priori 
uncertainty. 

Whan  a  sequence  of  measurements  are  taken  over  a  period  of 
time,  instantaneously  smallest  position  errors  may  not  give  the  smallest 
time  average  position  errors.  There  is  a  position-velocity  error  trade¬ 
off  that  relates  to  the  error  growth  between  measurements  as  well 
the  error  reduction  at  a  measurement.  This  question  is  dealt  with 
in  detail  for  a  simplified  model,  and  a  proposed  approach  is  given 
for  a  more  realistic  model. 

When  one  or  more  satellite  channels  are  lost  or  heavily 
jammed,  the  recursive  Kalman  filter  can  carry  on  with  a  smaller 
measurement  vector.  This  will  degrade  system  accuracy,  so  it  may 
be  advantageous  to  consider  substituting  backup  sensors  in  place 
of  the  missing  or  degraded  satellite  signals.  A  suboptima l  parti¬ 
tioned  state  vector  approach  to  the  problem  has  been  developed, 
along  with  soma  guidance  regarding  which  substitute  sensor  ought 
to  be  selected,  and  under  what  conditions  the  substitution  ought 
to  be  made.  The  answers  to  these  questions  depend  upon  the 
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cMWlMlon.  tMt  «l.t  —a.  eo»pon.nt.,  .cd  upoa  the 
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sensors  and  tfti  «Wlw«  'wMr*'  '  "  '  . 
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.  ,  the  optimal  static  measurement 

,vj.  v»ri*tv  can  be  imagined,  tne  p 

■election  problam,  the  dynamic  or  sequential  measurement  selection 

problem,  and  the  sensor-substitution  problem.  Suggested  direction 

for  further  study  and  simulation  tests  have  been  given 

these  areas. 
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